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The  small  cross  section  of  the  bridges  resulted  in  a  uniform  current 
distribution  In  the  bridges  at  all  temperatures,  thus  enabling  the  experi¬ 
mental  determination  of  the  critical  currant  versus  temperature  relation 
for  all  temperatures  between  1  j and  0.05  T.  where  ”T\ is  the  transition 
temperature  of  the  electrodes.  This  experimental  I  XT)  relation  was  found 
to  agree  with  theoretical  predictions  based  on  Usad£i's  equations.  This 
successful  fit  of  theory  to  experiment  produced  a  reduced  length  i  for 
each  sample;  the  samples  had  reduced  lengths  ranging  from  8.2  to  30. 

^>The  low-voltage  current-voltage  (I-V)  characteristics  were  studied  in 
detail.  It  was  found  that  these  I-V  curves  could  be  characterized  by  a 
temperiture  independent  relaxation  time  T  ^ .;  this  relaxation  time  was 
observed  to  be  proportional  to  £  .  At  low1 temperatures  the  I-V  curves 
were  hysteretlc,  with  the  curves  returning  to  the  zero-voltage  axis  at  a 
return  current  I_.  It  was  observed  that  1^  was  independent  of  temperature 
at  low  temperatures,  having  the  low  temperature  value  of. X^R  ~  5fi/2ei:  ... 
This  hysteresis  could  be  definitely  attributed  to  the  effective  relaxation 
time,  in  contrast  to  the  heating  or  capacitance  explanations  used  for 
all-superconductor  microbridges. 


At  all  temperatures  the  low-voltage  I-V  curves  were  straight  lines  for 
voltages  greater  than  S&/2eK^£. 

It  was  found  that  the  nature  of  th  I-V  curve  with  applied  microwave 
radiation  was  also  dependent  on  T  ...  If  the  microwave  period  was  shorter 
than  T  , .  or  li/2el  R,  there  was  aemcrowave-enhanced  critical  current.  For 
microwave  periods  longer  than  both  these  times  there  was  no  enhancement. 
This  result  indicated  that  the  critical-current  enhancement  was  most  likely 
caused  by  a  microwave-induced  spatial  broadening  of  the  quasiparticle  dis¬ 
tribution  in  the  bridge.  . 


Constant-voltage  steps  were  observed  in:.the  I-V  curves  for  all  frequencies 


of  microwave  radiation  used.  For  the  shortest  samples  steps  were  observed 


7 

-m 


at  voltages  greater  tban2'mV.  Subharmonic  steps  were  observed  whenever 
the  critical  current  was  enhanced. 


-^>The  high  resistance  of  these  samples  permitted  measurement  of  the  I-V 
curves  for  voltages  much  greater  than  the  energy  gap  of 'the  electrodes. 

A  high-voltage  peak  in  the  dV/dl  curves  was  observed,  and  it  was  speculated 
that  this  peak  was  due  to  an  extra  dissipation  in  the  electrodes.  In 
addition,  the  low-voltage  resistance  was  found  to  have  a  temperature  de¬ 
pendence  which  agreed  with  the  phenomenological  theory  of  Plppard  et  al. , 
a  theory  which  attributes  the  extra  resistance  to  energy  dissipation  in 

the  electrodes. 
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PROPERTIES  OF  SUPERCONDUCTOR-NORMAL  METAL- 
SUPERCONDUCTOR  MICROBRIDGES 

John  Michael  Warlaumont,  Ph.D. 

Cornell  University,  198C 

This  thesis  describes  the  fabrication  and  electrical  measurement  of 
short,  high  resistance  (R  ^  1  ohm)  superconductor-normal  metal - 
superconductor  microbridges. 

These  small  (length  and  width  as  small  as  200  nm) ,  planar  structures 
were  produced  using  electron-beam  lithography  and  ion-beam  etching.  Lead- 
copper-lead  structures  were  primarily  studied;  other  materials  were  also 
briefly  investigated. 

The  small  cross  section  of  the  bridges  resulted  in  a  uniform  current 
distribution  in  the  bridges  at  all  temperatures ,  thus  enabling  the 
experimental  determination  of  the  critical  current  versus  temperature 
relation  for  all  temperatures  between  T  and  0.05  T  ,  where  T  is  the 

C  v# 

transition  temperature  of  the  electrodes.  This  experimental  IC(T)  relation 
was  found  to  agree  with  theoretical  predictions  based  on  Usadel's 
equations.  This  successful  fit  of  theory  to  experiment  produced  a  reduced 
length  t  for  each  sample;  the  sampler,  had  reduced  lengths  ranging  from 
8.2  to  30. 

The  low-voltage  current-voltage  (I-V)  characteristics  were  studied 
in  detail.  It  was  found  that  these  I-V  curves  could  be  characterized  by  a 

temperature  independent  relaxation  time  t  this  relaxation  time  was 

2 

observed  to  be  proportional  to  L  .  At  low  temperatures  the  I-V  curves 
were  hysteretic,  with  the  curves  returning  to  the  zero-voltage  axis  at  a 
return  current  IR.  It  was  observed  that  IR  was  Independent  of  temperature 


at  low  temperatures,  having  the  low  temperature  value  of 
I^R  ^  5ti/2ereff.  This  hysteresis  could  be  definitely  attributed  to  the 
effective  relaxation  time,  in  contrast  to  the  heating  or  capacitance 
explanations  used  for  all -superconductor  microbridges. 

At  all  temperatures  the  low- voltage  I-V  curves  were  straight  lines 
for  voltages  greater  than  51f/2eref^. 

It  was  found  that  the  nature  of  the  I-V  curve  with  applied  micro- 
wave  radiation  was  also  dependent  on  t^.  If  the  microwave  period  was 
shorter  than  or*fi/2eIcR,  there  was  a  microwave-enhanced  critical 
current.  For  microwave  periods  longer  than  both  these  times  there  was 
no  enhancement.  This  result  indicated  that  the  critical-current 
enhancement  was  most  likely  caused  by  a  microwave-induced  spatial 
broadening  of  the  quasiparticle  distribution  in  the  bridge. 

Constant-voltage  steps  were  observed  in  the  I-V  curves  for  all 
frequencies  of  microwave  radiation  used.  For  the  shortest  samples  steps 
were  observed  at  voltages  greater  than  2  mV.  Subharmonic  steps  were 
observed  whenever  the  critical  current  was  enhanced. 

The  high  resistance  of  these  samples  permitted  measurement  of  the 
I-V  curves  for  voltages  much  greater  than  the  energy  gap  of  the  electrodes. 
A  high-voltage  peak  in  the  dV/dl  curves  was  observed,  and  it  was 
speculated  that  this  peak  was  due  to  an  extra  dissipation  in  the 
electrodes.  In  addition,  the  low-voltage  resistance  was  found  to  have  a 
temperature  dependence  which  agreed  with  the  phenomenological  theory  of 
Pippard  et  al_. ,  a  theory  which  attributes  the  extra  resistance  to  energy 
dissipation  in  the  electrodes. 


Chapter  1 
INTRODUCTION 

The  subject  of  this  work  is  the  experimental  investigation  of  a 
specific  kind  of  superconducting  weak  link— planar  superconductor- 
normal  metal -superconductor  microbridges.  But  the  goal  of  this  work  1 
not  only  to  characterize  the  physical  properties  of  SNS  microbridges, 
but  also  to  gain  a  better  understanding  of  the  other  types  of  super¬ 
conducting  junctions  by  noting  both  the  similarities  and  differences 
between  the  samples  of  this  work  and  other  systems.  In  fact,  we  have 
also  discovered  phenomena  which  are  applicable  to  physical  systems 
other  than  superconducting  junctions. 

In  this  introduction  a  very  brief  discussion  of  superconducting 
junctions  will  be  given,  as  well  as  an  explanation  of  how  this  work 
relates  to  the  previous  work  in  the  field.  The  remainder  of  this 
thesis  will  describe  in  detail  our  experimental  observations  and  our 
understanding  of  the  physical  processes  behind  these  observations. 

The  term  "superconducting  junction"  applies  to  any  conducting 
link  between  two  large  superconducting  electrodes;  the  electrodes  are 
assumed  to  behave  essentially  like  bulk  superconductors,  and  the  link 
between  them  is  assumed  to  be  sufficiently  weak  that  the  electro¬ 
magnetic  properties  of  the  sample  are  determined  mostly  by  the  link. 
The  basic  ideas  behind  such  structures  were  first  developed  by 
Josephson^  in  1952.  Since  that  pioneering  work  there  have  been 
thousands  of  papers  published  on  this  subject.  (A  recent  review  has 
been  given  by  Likharev.  )  But  in  spite  of  this  tremendous  amount  of 


effort,  there  are  still  a  large  number  of  unexplained  experiments, 
particularly  those  which  involve  non-equilibrium  situations  or  time- 
dependent  phenomena.  The  major  difficulty  on  the  theoretical  side  has 
been  the  very  complicated  microscopic  theory  of  superconductivity.  On 
the  experimental  side,  it  has  been  very  difficult  to  fabricate  well- 
defined  junctions  due  to  the  very  small  lengths  involved  and  the  often- 
undesirable  mechanical  and  chemical  properties  of  superconducting 
materials.  We  have  attempted  to  simplify  the  theoretical  analysis  by 
making  samples  which  have  a  very  well-defined  weak  region;  this  has  been 
accomplished  by  using  a  normal  metal  as  the  bridge  between  the  super¬ 
conducting  electrodes.  In  addition,  we  have  taken  advantage  of  modern 
lithographic  technology  to  significantly  reduce  the  size  of  the  bridge. 
We  have  succeeded  in  producing  samples  which  have  been  of  sufficient 
quality  to  allow  verification  of  some  previously  unproved  theories;  in 
addition,  we  have  observed  phenomena  which  have  not  been  anticipated  by 
any  of  the  available  theories. 

There  are  several  different  types  of  superconducting  junctions,  and 

only  one  type,  the  thin-film  microbridge,  will  be  discussed  here.  The 

3  4 

reader  is  referred  to  the  works  of  Jackel  and  Henkels  for  descriptions 
of  the  other  types. 

A  schematic  diagram  of  a  thin-film  superconducting  microbridge  is 
given  in  Figure  1.1.  The  structure  consists  of  two  superconducting 
electrodes  (banks)  connected  by  a  small  metal  link  (bridge).  The  metal 
films  are  deposited  on  an  insulating  substrate,  and  they  are  sculptured 
into  the  desired  shape  by  a  variety  of  techniques.  They  can  be  scribed 


L  -  3000  A 
w  ~  3000  A 
3b  1500  A  thick 
:u  —  600  A  thick 


Figure  1.1  Schematic  diagram  of  an  SNS  microbridge. 
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with  a  metal  stylus  ,  patterned  by  depositing  the  films  through  a 
mechanical  mask6,  or  they  can  be  tailored  using  lithographic  techniques. 
The  latter  method  was  used  almost  exclusively  in  this  work.  The  type  of 
microbridge  can  be  classified  further:  if  the  bridge  region  is  composed 
of  a  normal  metal,  it  is  called  a  superconductor-normal -superconductor 
(SNS)  microbridge;  if  the  bridge  is  thinner  than  the  banks,  it  is  a 
variable-thicKness  microbridge  (VTB);  if  both  the  bridge  and  the  banks 
are  the  same  thickness  and  the  same  material,  it  is  a  constriction-type 
(ScS)  microbridge;  and  if  the  weak  region  is  formed  by  underlying  the 
weak  region  with  a  lower  T  material,  it  is  a  proximity-effect  micro- 
bridge.  There  is  one  important  point  to  be  made  here:  SNS  microbridges 
and  proximity-effect  microbridges  are  not  the  same.  In  the  proximity- 
effect  structure  the  bridge  has  a  finite  transition  temperature,  usually 
above  IK.  In  the  SNS  microbridge,  the  bridge  does  not  become  super¬ 
conducting  at  any  temperature. 

1.1  Basic  Theory  of  Microbridges 

All  of  the  basic  phenomena  of  microbridges  are  the  result  of  one 
fact:  the  order  parameters  of  the  opposed  electrodes  are  coupled 
together  in  amplitude  and  in  phase  by  the  weak  link.  It  is  the  nature  of 
this  coupling  which  determines  the  properties  of  any  given  junction. 

1.1.1  Critical -current  and  current-phase  relation 
If  the  current  through  the  junction  is  less  than  a  critical  super¬ 
current  Ic,  the  coupling  between  the  electrodes  can  be  completely 
described  by  the  current-phase  relation  f{<j>).  In  this  case  the  current 
is  carried  compl etely  by  the  supercurrent  I  ,  where  I  is  given  by 
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I,  *  lcfM- 

For  an  Ideal  junction,  f(<t>)  «  sin$.  For  a  non-ideal  junction,  f(<j>) 
can  deviate  from  the  sinusoidal  form,  but  it  always  has  two  properties: 
f(4>)  is  always  an  odd  function,  and  f (4>)  is  always  periodic  with  period 
M  *  2ir;  f{$)  is  not  necessarily  single-valued.  Jackel  has  made  exten¬ 
sive  measurements  of  f <4>)  for  several  different  types  of  microbridges, 
and  his  work  contains  an  excellent  discussion  of  the  results. 

The  current-phase  relation  (CPR)  of  SNS  microbridges  has  been 
theoretically  studied  in  great  detail.  These  theories  are  described  in 
Chapter  4;  they  all  have  essentially  the  same  result.  Except  for  very 
short  bridges  at  very  low  temperatures,  the  CPR  should  be  sinusoidal. 

The  microbridges  of  this  work  all  had  measured  CPRs  which  were  sinusoidal, 
in  agreement  with  the  theories.  Thus,  as  far  as  the  zero-voltage  CPR  is 
concerned,  our  microbridges  were  ideal. 

The  value  of  I  depends,  of  course,  on  the  cross-sectional  area  of 
the  weak  link.  For  sufficiently  small  links,  I£  is  inversely  proportional 
to  the  normal-state  resistance  R  of  the  bridge;  thus,  to  eliminate  the 
areal  dependence  of  Ic,  the  value  of  ICR  is  usually  reported.  For  an 

ideal  microbridge,  with  banks  and  bridge  made  from  the  same  material,  the 

2 

I£R  product  close  to  Tc  should  be  given  by 

I  R  *  a(T  -T)  where  a  £  635  yV/K. 

w  C 

At  low  temperatures,  the  I  R  product  should  be  independent  of  tempera- 

v 

ture,  attaining,  for  small  microbridges  in  the  dirty  limit,  a  maximum  at 

2 

T  *  0  of  316  T  microvolts,  where  T  is  in  degrees  K. 

C  w 

For  an  SNS  structure,  on  the  other  hand,  close  to  T  the  I  R  product 

v  w 
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should  be  proportional  to  (T  -T)  ;  at  low  temperatures  I  R  should  be 

c  c 

nearly  exponential  with  temperature.17  We  have  compared  the  temperature 
dependence  of  I  R  for  our  samples  to  the  predictions  of  the  theories 

w 

applicable  to  these  samples,  and  we  have  found  very  good  agreement.  In 
fact,  for  T  <  0.6  Tc,  our  results  agree  both  in  form  and  in  magnitude  with 
the  theory. 

The  above  temperature  dependences  of  I£R  are  valid  only  for 

sufficiently  small  microbridges.  For  all -superconductor  structures,  if  the 

bridge  becomes  much  longer  than  the  coherence  length,  the  temperature 

2 

dependence  changes  radically  ;  this  is  not  the  case  for  SNS  structures. 

If  the  width  of  an  all -superconductor  microbridge  is  larger  than  the 
penetration  depth,  the  ICR  product  will  be  reduced  due  to  the  formation  of 
vortices  inside  the  bridge.  This  can  also  happen  in  SNS  junctions;  in 
fact,  nearly  all  of  the  previous  work  on  SNS  layered  junctions  produced 
ICR  products  which  were  limited  at  low  temperatures  due  to  the  relatively 
large  cross-sectional  area  of  the  bridge.  This  was  not  a  problem  in  the 
SNS  structures  of  this  work,  since  for  our  samples  the  penetration  depth 
is  of  the  order  of  1  mm7,  which  is  much  larger  than  the  width  or  thick¬ 
ness  of  the  bridges  used  here. 

1.1.2  Low- voltage  phenomena 

If  a  voltage  V  is  present  across  the  bridge,  there  are  two  addi¬ 
tional  factors  which  must  be  considered.  First  of  all,  in  addition  to 
the  supercurrent  there  is  a  normal  current  through  the  junction,  given 
by  In  *  V/R.  Second,  the  phase  difference  between  the  banks  is  now  time- 
varying,  with  the  variation  given  by  d$>/dt  s  (2e/-K)V.  Thus,  for  an  ideal 
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junction,  the  current  through  the  junction  is  given  by 

I  *  Is  +  In  *  Ics1n<J>  +  V/R. 

This  model  of  an  ideal  junction  is  often  called  the  res istively -shunted 
junction  (RSJ)  model.  For  this  model  the  current  versus  voltage  ( I -V ) 
characteristic  is  given  by  Vdc  *  R(Idc»I^)1//2. 

The  problem  with  the  RSJ  model  is  that  no  microbridge  studied  has 
ever  had  the  I-V  curve  predicted  by  the  RSJ  model.  The  I-Vs  of  real 
microbridges  often  show  hysteresis,  bumps,  curvature  in  the  wrong 
direction,  and  other  miscellaneous  features.  Most  of  these  deviations 

from  the  RSJ  model  are  caused  by  heating  or  by  relaxation  time  effects. 

2 

These  effects  have  been  described  in  detail  by  Likharev. 

It  will  be  shown  in  Chapter  4  that  our  SNS  junctions  did  not  follow 
the  predictions  of  the  RSJ  model.  It  will  be  shown  that  joule  heating 
in  the  bridge  was  not  important.  Instead,  it  was  found  that  these  SNS 
structures  had  an  intrinsic  relaxation  time  which  caused  significant 
hysteresis  in  the  I-V.  We  observed  that  the  hysteresis  could  be 
characterized  by  the  IR  product  at  which  the  I-V  returned  to  the  zero- 
voltage  axis.  This  characteristic  IR  product  was  independent  of  tempera¬ 
ture  at  low  temperatures,  and  the  low- temperature  value  was  propor¬ 
tional  to  the  square  of  the  reduced  length  of  the  bridge.  Both  of 
these  phenomena  are  predicted  by  time-dependent  Ginzburg-Landau 
calculations  for  SNS  structures. 

In  addition,  we  found  that  the  observed  low-voltage  curvature  in 
the  dc  I-V  could  be  correlated  with  the  intrinsic  relaxation  time, 
although  precise  comparisons  with  theory  could  not  be  made. 
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1.1.3  RF  response  of  microbridges 

The  fact  that  the  supercurrent  oscillates  in  time  when  there  is  a 
voltage  across  the  bridge  has  led  to  a  great  number  of  experiments  in 
which  the  microbridge  is  driven  by  an  ac  current  as  well  as  a  dc  current. 
Such  experiments  have  produced  two  results.  The  first  is  that  for  small 
ac  currents  the  critical  current  is  increased  above  the  value  it  had  with¬ 
out  an  ac  drive  current.  This  enhancement  has  been  studied  intensively 
in  the  last  few  years,  and  our  SNS  results  have  led  to  a  few  surprises. 
For  all  superconductor  microbridges  it  is  widely  believed  that  the  ac 
field  enhances  the  energy  gap  in  the  bridge,  thereby  increasing  the 
critical  current.  However,  the  results  of  this  work  indicate  that  the 
enhancement  in  SNS  structures  is  due  to  a  spatial  redistribution  of  the 
electrons  in  the  bridge,  rather  than  an  increase  in  an  energy  gap.  The 
evidence  for  this  explanation  lies  in  the  very  strong  correlation  between 
the  effective  relaxation  time  and  the  frequency  of  ac  current  required 
to  increase  the  critical  current. 

The  second  phenomena  observed  when  an  ac  current  is  applied  is  the 
observance  of  constant-voltage  steps  in  the  I-V  curve.  These  steps  occur 
at  voltages  which  are  multiples  of  the  characteristic  voltage  (ii/2e)tJrf . 
(Henkels*  has  given  a  good  review  of  this  phenomena.)  In  this  work,  as 
well  as  in  all  previous  works  on  microbridges,  steps  always  occurred  at 
multiples  of  this  characteristic  voltage.  But  for  some  frequencies 
steps  also  appeared  at  submultiples  of  this  characteristic  voltage.  The 
minimum  frequency  at  which  these  subharmonic  steps  occurred  in  this  work 
was  proportional  to  the  effective  relaxation  time  determined  from  the 
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reduced  length  of  the  bridge. 

The  major  disagreement  with  the  theory  of  rf-induced  steps  was  in  the 
magnitude  of  the  steps.  For  frequencies  above  a  frequency  proportional  to 
the  effective  relaxation  time,  there  was  strong  disagreement. 

In  summary,  we  have  discovered  that  the  effective  relaxation  time 
determined  from  dc  measurements  was  very  important  in  predicting  the  ac 
response  of  our  microbridges. 

1.1.4  High-voltage  effects 

The  relatively  high  resistance  of  our  samples  allowed  us  to  measure 
the  I-V  characteristic  at  voltages  above  the  gap  voltage  of  the  super¬ 
conducting  electrodes;  this  regime  had  not  been  reached  in  previous 
studies  of  SNS  structures.  We  found  that,  in  contrast  to  the  low-voltage 
measurements,  the  effect  of  the  normal  metal  bridge  on  the  electrodes 
was  also  important,  in  addition  to  the  effect  of  the  banks  on  the  bridge. 
It  was  observed  that  the  dV/dl  curve  had  features  which  could  be 
attributed  to  energy  dissipation  in  the  banks.  Previous  experimental 
workers  have  observed  high-voltage  structure  in  the  I-V  curves,  but  no 
previous  experimental  or  theoretical  work  has  considered  the  possibility 
of  an  extra  high-voltage  dissipation  in  the  banks. 

1 .2  Sample  Fabrication  and  Measurement 

The  very  small  length  scale  of  superconductivity  requires  samples 
to  be  very  small.  To  be  more  precise,  the  coherence  length  for  a  normal 
metal  is  given  by  £C(T)  =  (tfVfA/6‘rckT)  '  where  A  is  the  electron  mean 
free  path.  For  the  copper  films  used  in  this  work,  A  ~  15  nm.  (A  was 
not  measured  directly;  it  was  deduced  by  fitting  the  ICR  versus  T  data 
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to  the  theory.  A  direct  calculation  from  the  resistance  R  and  the 
physical  dimensions  of  sample  Pb-Cu-13  (R  =  .3  ohms,  length  L  %  200  nm, 
width  w  S'  200  nm,  thickness  d  S'  60  nm)  gave  a  value  of  A  =  25  nm.  But 
this  direct  calculation  is  only  approximate  due  to  the  uncertainty  in 
the  actual  dimensions  of  the  bridge.)  This  means  that  for  T  =  7K  (the 
transition  temperature  of  the  lead  banks)  £0(T=7K)  S'  40  nm.  Now,  in  order 
to  get  useful  samples  we  needed  bridges  whose  lengths  were  of  the  order 
of  10  £0(T=7K)  or  smaller.  To  achieve  this,  we  could  not  use  the  tra¬ 
ditional  photo-lithography;  instead,  we  developed  our  own  version  of 
electron-beam  lithography.  With  this  technique,  we  were  able  to  achieve 
bridges  of  length  L  S' 8  $0(T=7K). 

The  fact  that  we  were  interested  in  samples  composed  of  different 
materials  led  to  difficulties  in  obtaining  a  good  metallic  contact  between 
the  different  materials.  This  necessitated  the  introduction  of  a 
cleaning  step  (ion-bombardment)  between  the  separate  metal  depositions. 

We  were  also  interested  in  various  materials:  tin,  lead,  niobium, 
aluminum,  indium,  and  copper.  To  this  end,  both  sputtering  and  thermal 
evaporation  were  used  to  deposit  the  films . 

All  together,  several  different  combinations  of  materials  were  tried. 
The  best  results  were  obtained  using  lead  for  the  superconducting  elec¬ 
trodes,  and  copper  for  the  normal -metal  bridges.  Only  the  results  of  these 
Pb-Cu-Pb  structures  will  be  presented  here  in  detail.  The  results  on 
the  other  structures  will  be  briefly  presented  in  Chapter  5. 

The  measurements  performed  were  quite  standard,  at  least  in  our 
laboratory.  A  discussion  of  the  measurement  techniques  will  therefore  be 
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brief,  with  only  those  techniques  novel  to  this  work  presented  in 
detail . 
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Chapter  2 
SAMPLE  FABRICATION 

One  of  the  most  time-consuming  parts  of  this  work  was  the  preparation 
of  samples.  Thin  film  samples  were  required  in  which  the  planar  dimen¬ 
sions  were  very  small,  of  the  order  of  200  nm.  Structures  this  small 
could  not  be  reliably  produced  using  photolithography,  so  an  electron- 
beam  lithography  (EBL)  system  was  developed;  this  system  could  produce 
structures  with  planar  dimensions  as  small  as  150  nm.  In  addition, 
several  processing  steps  often  occurred  between  film  depositions,  so  that 
ion-beam  etching  or  sputter  etching  had  to  be  employed  between  depositions 
in  order  to  achieve  good  metallic  contact  between  successively  deposited 
films. 

The  sequence  of  fabrication  steps  was  usually  the  following:  a 
gross-electrode  structure  was  put  onto  a  cleaned  substrate  using  photo¬ 
lithography.  The  bridge  component  of  the  microbridge  was  then  deposited 
using  EBL  and  the  lift-off  method.  Next,  the  banks  component  was  laid 
down,  also  with  EBL  and  the  lift-off  method.  (In  the  case  of  niobium 
banks,  the  banks'  geometry  was  sculptured  with  sputter  etching.)  Finally, 
the  gross  electrode  structure  was  expanded  and/or  thickened,  and  the 
sample  was  ready  for  testing.  This  chapter  is  devoted  to  a  detailed 
discussion  of  the  various  steps  of  this  process. 

2.1  Substrate-Selection  and  Cleaning 

The  first  substrates  employed  were  thin  (.006  in)  glass  cover 
slips.  These  substrates  did  not  prove  to  be  satisfactory  for  two 
reasons:  they  were  so  thin  that  they  broke  very  easily,  and  glass  has 
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a  very  poor  thermal  conductivity.  Good  thermal  conductivity  is  an 
important  feature  in  a  substrate  since  it  aids  in  the  cooling  of  the 
substrate  during  film  depositions  and  it  improves  the  performance  of 
the  microbridge  laid  upon  it,  although  no  experiments  were  done  to 
meausre  the  importance  of  this  latter  advantage.  Fused  quartz  was  also 
tried  (briefly)  as  a  substrate,  but  although  it  has  adequate  strength 
it  has  a  poor  thermal  conductivity.  The  final  choice  was  sapphire,  and 
it  worked  quite  well  as  a  substrate  material.  Single  crystal  sapphire 
is  very  strong,  and  its  thermal  conductivity  is  of  the  same  order  of 
magnitude  as  that  of  a  metal  at  all  temperatures  of  interest. 

Single  crystal  sapphire  substrates  were  obtained  from  Adolf  Meller 
Co.  with  a  standard  polish  on  one  side  and  an  extra-fine  polish  on  the 
side  used  for  the  sample.  They  came  pre-cleaned,  but  an  ultrasonic 
cleaning  in  isopropyl  alcohol  was  necessary  to  remove  any  dust  on  the 
surface.  Generally  each  substrate  could  be  used  only  once,  since  some 
of  the  processing  steps  left  resist  residue  which  was  impossible  to 
remove  completely  using  any  available  cleaning  technique. 

2.2  Photolithography  and  the  Lift-off  Technique 

Most  of  the  gross  electrode  structures  were  obtained  using  photo¬ 
lithography  and  the  lift-off  technique.  This  method  has  been  in  use 

4 

for  many  years  and  a  good  explanation  of  it  has  been  given  by  Henkel s, 
but  a  brief  discussion  of  the  method  is  useful  here.  (Figure  2.1 
depicts  the  method.) 

First,  a  thin  (200  nm  to  1000  nm)  film  of  photo  resist  was  spun 
onto  the  substrate.  (The  substrate  could  have  been  clean,  or  it  could 
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have  had  previously-defined  thin  film  structures  on  it,)  Shipley  AZ1350B 
and  Shipley  AZ13500  positive  resists  were  used.  After  application,  the 
resist  was  baked  at  80°C  for  ten  minutes.  A  mask  (either  a  film  negative 
or  a  metal  film  on  a  transparent  backing)  was  laid  over  the  top,  and  the 
mask  was  held  in  intimate  contact  with  the  substrate  in  a  vacuum  chuck. 
Ultraviolet  light  was  then  shone  through  the  mask.  A  high  pressure 
mercury  lamp  with  an  incident  flux  of  about  250  foot-candles  gave 
exposure  times  of  about  three  minutes.  After  the  exposure,  the  mask  was 
removed  and  the  resist  was  developed  in  an  aqueous  developer;  this 
developer  removed  the  exposed  areas  of  the  resist. 

After  the  resist  stencil  was  completed,  the  lift-off  method  was 
employed:  the  desired  metal  film  was  evaporated  over  the  top  of  the 
resist,  then  the  resist  was  removed  using  a  resist  stripper  (acetone); 
metal  which  was  on  top  of  the  resist  was  removed  along  with  the  resist.  ■■ 
Remaining  on  the  substrate  was  the  desired  result— a  metal  film  whose 
geometry  was  the  same  as  that  of  the  transparent  area  of  the  mask. 

Although  the  photolithographic  and  lift-off  techniques  were 
relatively  simple,  they  had  significant  drawbacks.  The  wavelength  of 
the  ultraviolet  light  limited  the  resolution  to  about  500  nm,  a 
resolution  which  was  not  fine  enough  for  the  desired  microbridges; 
electron-beam  lithography  was  used  to  obtain  finer  resolution.  In 
addition,  the  lift-off  method  was  incompatible  with  the  sputtering 
process  used  for  the  deposition  of  niobium  films;  in  these  cases  it  was 
necessary  to  use  a  sputter-etching  technique  to  define  the  niobium 
films. 


2.3  Electron-Beam  Lithograph: 


Q 

Electron-beam  lithography  (EBL)  is  quite  similar  to  photolithography; 
only  the  method  of  exposure  is  substantially  different.  In  EBL,  the 
resist  is  exposed  by  electrons  in  a  scanning  electron  microscope  (SEM). 

The  resolution  of  EBL  is  determined  not  by  the  wavelength  of  the  electrons 
(the  electrons  used  have  a  wavelength  of  about  .01  nm)  but  by  the  profile 
of  the  electron  beam  in  the  resist.  The  factors  which  determine  this 
profile  will  be  discussed  later.  The  minimum  resolution  achieved  in  this 
work  was  about  150  nm. 

2.3.1  Electron-beam  resist 

The  resist  used  in  this  work  was  Dupont  Elvacite  2041  (PMMA),  a 
highly  crossed-1 inked  polymer  (polymethylmethacrylate)  with  a  mean 
molecular  weight  of  about  500,000.  The  properties  of  this  resist  have 
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been  described  in  detail  by  Wolf  al_.  The  resist  was  obtained  in 
powder  form;  the  powder  was  dissolved  in  methyl  isobutyl  ketone  (MI K) , 
filtered  through  a  Mi 11 i pore  filter,  and  spun  onto  a  substrate  at  about 
5000  rpm.  It  was  found  that  a  4-1/25?  (by  weight)  concentration  of  PMMA 
in  MIK  gave  a  resist  thickness  of  about  300  nm;  an  8%  solution  gave  a 
thickness  of  about  600  nm.  After  application,  the  resist  was  baked  at 
100°C  for  one  hour.  (Baking  at  160°C  for  one-half  hour  gave  better 
adhesion  and  slightly  better  exposure  characteristics,  but  this  baking 
temperature  proved  to  be  damaging  to  metal  films  which  may  have  been 
present  under  the  resist.)  The  8%  solution  proved  to  be  markedly  non- 
uniform,  so  thicker  resist  coatings  were  produced  by  applying  more  than 
one  coat  of  the  4-1/2%  solution;  the  samples  were  baked  for  one-half  hour 
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between  coats  and  for  one  hour  after  the  final  coat.  Two  coats  of  4-1/2% 
solution  gave  a  thickness  of  about  550  nm  and  three  coats  of  4-1/2%  solution 
gave  a  thickness  of  about  800  nm.  After  baking,  10  nm  of  aluminum  was 
evaporated  onto  the  sample  to  prevent  charging  in  the  SEM. 

After  exposing  the  resist  in  the  SEM,  the  resist  was  developed  with 
the  following  procedure: 

Soak  in  0.1  Molar  NaOH  to  remove  A1  coating. 

Rinse  in  distilled  water. 

Develop  in  3:1  isopropyl  alcohol :MIK  for  30  sec  to  5  minutes. 

Ultrasonic  3:1  isopropyl  alcohol :MIK  for  5  sec. 

Rinse  in  isopropyl  alcohol. 

Spin  dry. 

Because  PMMA  is  a  positive  resist,  the  areas  of  the  resist  exposed 
to  the  electron  beam  were  removed  by  the  development  process.  After 
development,  the  samples  were  ready  for  the  deposition  step. 

2.3.2  Electron-beam  exposure  system 

Resolution  _ _ 

The  electron-beam  resist  coatings  were  exposed  in  an  AMR  Model  1000 
scanning  electron  microscope  (SEM).  This  microscope  had  a  minimum  probe 
diameter  of  about  20  nm,  but  its  usual  resolution  was  about  50  nm.  This 
resolution  was  sufficient  for  all  the  samples  produced,  since  the 
minimum  pattern  size  in  the  resist  was  governed  by  the  scattering  of 
electrons  in  the  resist.  The  electrons  had  an  energy  of  20  keV,  and  at 
this  energy  the  scattering  in  the  resist,  neglecting  substrate  effects, 
limited  the  resolution  of  the  exposed  pattern  to  about  one-fourth  the 
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thickness  of  the  resist.  If  the  substrate  had  been  coated  with  a  moderate 

to  high  density  material,  this  minimum  feature  size  could  be  increased  to 

about  one-half  the  resist  thickness  (or  larger).  ®  This  scattering  limit 

could  have  been  reduced  by  using  an  SEM  with  a  higher  accelerating 
g 

voltage. 

Registration 

PMMA  has  a  relatively  low  sensitivity  to  electrons;  as  a  result,  the 
EBL  system  could  only  expose  relatively  small  areas.  (An  area  1  mm  by 
1  mm  square  required  a  minimum  of  one  hour  to  expose.)  Much  larger  areas 
were  required  for  making  the  electrical  contacts  necessary  for  testing. 

For  some  of  the  samples,  the  gross  electrode  structure  was  put  on  the  sub¬ 
strate  orior  to  the  EBL  steps.  In  other  cases,  ones  in  which  the  necessary 
electrodes  could  not  survive  the  EBL  procedure,  small  silver  pads  were  put 
on  the  substrate  to  aid  in  the  registration  of  the  EBL  exposure;  silver 
was  used  because  it  was  easy  to  evaporate,  stood  up  to  the  EBL  processing 
steps,  and  has  a  high  electron  cross  section,  thus  providing  a  good  image 
in  the  SEM.  These  silver  pads  were  especially  useful  in  the  banks'  expos¬ 
ure,  since  the  thin  bridge  lines  were  very  hard  to  locate  without  these 
registration  pads.  Figure  2.2  shows  a  few  of  the  gross  electrode/ 
registration  structures  which  were  used. 

These  registration  markers  were  also  very  useful  for  focusing.  The 
usual  procedure  was  to  position  the  sample  at  low  magnification,  move 
slightly  away  with  either  the  electrostatic  beam  shift  or  the  mechanical 
stage  translation,  go  to  high  magnification  to  focus,  reposition  at 
low  magnification,  and  then  expose  at  high  magnification.  It  was  not 
possible  to  focus  directly  at  the  place  to  be  exposed,  since  this 
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Figure  2.2  Registration  pads  and  electrode  structures  for  three  samples 
Top  -  Nb-Al-4.3;  middle  -  Nb-Cu-1/4;  bottom  -  Pb-Cu~13. 

The  gross  electrodes  were  patterned  using  an  aluminum- foil 
stencil;  the  electrodes  and  registration  pads  were 
lithographically  produced. 
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focusing  would  expose  the  resist  in  an  unwanted  pattern.  The  focusing 
had  to  be  done  close  to  the  desired  exposure  location  since  the  SEM  could 
not  keep  its  resolution  over  a  large  area  of  the  sample.  If  there  wasn't 
a  registration  pad  close  enough  to  focus  on,  a  minute  or  so  of  exposure 
at  the  highest  magnification  (100,000  X}  would  blister  the  resist  and 
thus  create  a  defect  to  focus  on. 

Line  exposure 

The  SEM  had  a  fair  amount  of  noise  (mostly  at  60  Hz)  in  the  scan 
generator-deflection  coil  circuits.  Because  of  this,  the  finest  pattern 
which  could  be  drawn  was  a  line  in  the  resist,  exposed  in  a  single  pass. 

For  many  purposes  (such  as  the  bridge  line  exposures  or  the  gap  patterns 
for  sputter  etching)  a  single  line  was  quite  satisfactory.  Such  lines 
were  drawn  with  the  SEM  in  line-scan  mode.  The  beam  was  blanked  except  for 
a  portion  of  a  single  line  scan;  it  was  important  to  blank  both  the  begin¬ 
ning  and  end  of  a  line  in  order  to  avoid  the  irregularities  which  the 

scan  generator  put  at  either  end.  In  this  mode,  the  necessary  exposure 

•8  —1 0 

was  about  1.3x10  .  coul/cm.  (At  a  specimen  current  of  1.3  x  10"  J  amps, 

a  .02  cm  line  could  be  drawn  in  2  sec.) 

Rectangle  exposure 

The  micrcbridges  which  were  fabricated  were  very  small;  the  bridge 

region  which  required  the  fine  resolution  occupied  an  area  of  about  100 

square  microns.  Since  the  gross  electrode  structures  had  gaps  of  about 

2 

1  mm,  it  was  usually  necessary  to  expose  large  rectangles  (M  mm  )  which 
required  only  re’ativei y  poor  resolution.  This  was  most  easily  done  at 
low  magnification  in  the  partial  scan  mode  of  the  SEM.  Relatively  high 


21 


beam  currents  were  used  (4  x  lO”^  amps)  and  gaps  of  about  5  microns  could 
be  left  between  rectangles.  Figure  2.3  shows  the  use  of  such  rectangles 
in  the  production  of  a  Pb-Cu-Pb  microbridge.  The  beam  blank  was  not 
required  during  the  exposure  of  these  large  rectangles. 

Pattern  exposure— the  flying  spot  scanner 

Patterns  which  required  detail  finer  than  5  microns  and  which  could 
not  be  constructed  from  a  few  single  line  scans  were  exposed  with  the  SEM 
in  its  raster  mode.  In  this  mode  the  beam  scanned  the  sample  in  a 
raster  consisting  of  a  number  of  horizontal  lines.  The  pattern  was 
generated  by  blanking  the  beam  when  it  covered  areas  which  were  not  to 
be  exposed.  This  selective  blanking  of  the  beam  was  accomplished  with  a 
flying  spot  scanner  (FSS),  depicted  in  Figure  2.4. 

A  film  mask  was  made  in  the  same  geometry  as  the  desired  pattern, 
but  on  a  greatly  magnified  scale.  (The  magnification  of  the  mask  was 
the  SEM  magnification  at  which  the  exposure  was  to  be  done;  i.e.,  if  the 
exposure  was  to  be  done  at  10,000  X,  one  centimeter  on  the  mask  would  be 
one  micron  on  the  sample;  the  record  screen  had  dimensions  of  10  cm  by 
10  cm.)  This  mask  was  then  placed  over  the  record  screen  of  the  SEM,  and 
a  photomultiplier  was  used  to  detect  whether  or  not  the  record  screen 
could  be  seen  through  the  mask.  If  the  record  scope  trace  was  behind  a 
transparent  region  of  the  mask,  the  SEM  electron  beam  was  left  on;  if  the 
record  trace  was  behind  an  opaque  area  of  the  mask,  the  SEM  beam  was 
blanked.  Since  the  trace  of  the  record,  scope  was  slaved  to  the  same 
raster  as  the  SEM  beam,  the  SEM  beam  was  blanked  according  to  the  mask 
pattern,  and  the  resist  was  exposed  in  the  mask  pattern,  at  a  reduction 
corresponding  to  the  magnification  of  the  SEM  used  during  the  exposure. 
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Figure  2. 


3  Diagram  of  the  procedure  for  generating  electrode 

structures  using  electron-beam  lithography.  Figure  2.2 
(bottom)  shows  an  example  of  the  results  cf  this  proceoure. 
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For  example,  if  the  SEM  was  used  at  a  magnification  of  10,000  X,  the 
pattern  was  reduced  by  a  factor  of  10,000  in  a  single  step. 

The  circuit  used  to  interface  the  photomultiplier  to  the  beam  blank 
on  the  SEM  is  shown  in  Figures  2.5  and  2.6.  The  SEM  required  0  volts  to 
have  the  beam  on,  and  +5  volts  to  blank  the  beam.  This  circuit  contained 
a  threshold  detector  to  eliminate  the  zero  offset  and  noise  present  in  the 
photomultiplier  output.  It  also  contained  the  options  of  blanking  the 
beam  completely,  or  of  leaving  the  beam  completely  unblanked.  The  circuit 
was  also  designed  tc  duty  cycle  the  beam  while  it  was  in  the  "beam  on" 
state  in  order  to  reduce  the  current  in  the  beam;  this  feature  was  seldom 
used  since  it  had  adverse  effects  on  the  resolution.  The  photomultiplier 
was  run  with  about  one  kilovolt  between  the  first  dynode  and  the  anode, 
and  a  standard  resistor  network  was  used  for  biasing  the  intermediate 
dynodes.  The  best  performance  of  the  flying  spot  scanner  was  achieved 
with  the  record  scope  contrast  at  a  minimum. 

The  best  resolution  was  obtained  with  an  SEM  magnification  of 
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10,000  X.  At.  this  magnification,  a  sample  current  of  about  9x10  amps 
was  used.  The  horizontal  lines  had  a  scan  time  of  .013  sec/line;  the 
frame  scan  time  was  24  sec.  In  this  mode,  as  well  as  in  the  rectangle- 
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exposure  mode,  the  PMMA  had  a  sensitivity  of  about  1.8  x  10  C/cmw. 

The  most  important  consideration  in  the  design  of  the  flying  spot 
scanner  was  its  resolution;  the  resolution  was  characterized  by  the 
minimum  mask  dimension  which  could  be  resolved.  This  resolution  was 
determined  by  both  geometrical  and  response- time  factors.  There  were 
three  geometrical  factors: 


Duty  cycler 
output 


output 

Figure  2.5  Schematic  diagram  of  the  flying-spot  scanner  electronics. 


26 


Figure  2. 


Duty  Cycler  for  FSS 


n  n  .  -+  5V  on 

output  is  - 1  I - !  I — I  I - 0  off 

pulse  freq  *  0.06  or  0.6  MHz 

pulse  width  is  400  ns  to 

IC  is  Foirchiid  96L02C 


6  Schematic  diagram  of  the  duty-cycler  electronics  for 
the  FSS. 
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Resolution  of  the  film  mask— using  Kodalith  graphics  film  gave  a 
resolution  of  about  50  microns. 

Resolution  of  the  record  scope— the  SEM's  record  scope  had  a  resolu¬ 
tion  of  50  microns  or  less,  as  long  as  the  intensity  was  not  set  too 
high. 

Parallax  caused  by  the  fact  that  there  was  about  1/4  in  of  glass 
between  the  phosphor  on  the  inside  of  the  record  screen  and  the  film 
mask— placing  the  photomultiplier  24  in  from  the  mask  and  using  a  one 
inch  aperture  at  the  photomultiplier  reduced  this  error  to  a  few 
microns.  Parallax  also  caused  quite  a  bit  of  distortion  in  the  pattern, 
but  this  distortion  was  not  important  in  the  structures  used  in  this 
work.  The  use  of  lenses  between  the  mask  and  the  record  screen  and 
between  the  photomultiplier  and  the  mask  could  have  significantly  reduced 
the  distortion  as  well  as  eliminated  the  need  for  such  a  large  spacing 
between  the  photomultiplier  and  the  mask. 

These  geometrical  limits  were  quite  small;  the  significant  limita¬ 
tion  on  the  FSS  resolution  was  due  to  response-time  considerations.  A 
response  time  in  the  system  affected  the  resolution  because  if  the  system 
took  a  time  tf  to  respond  to  changes  in  light  intensity,  there  would  be 
a  resolution  limit  of  ds  -  t  dx/dt  where  dx/dt  is  the  speed  of  the 
record  scope  beam  across  the  record  screen.  There  were  two  response 
times  of  interest  in  the  system: 

Delays  in  the  photomultiplier,  in  the  threshold  detector  and  in 
the  SEM's  beam  blank  circuitry— the  total  of  these  delays  was  about  one 
microsecond.  This  delay  was  much  smaller  than  the  other  response  time, 
and  thus  posed  no  significant  limitation.  In  fact,  the  major  effect  of 
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this  delay  was  a  lateral  shift  in  the  pattern;  this  was  a  form  of  distor¬ 
tion  and  was  easily  corrected. 

Record  screen  phosphor  decay  time— this  was  the  significant  limita¬ 
tion  on  the  FSS  resolution.  The  phosphor  used  (type  Pll)  had  a  decay 
time  of  about  50  microseconds.  At  a  sample  current  of  10"^  amps  and  a 
magnification  of  10,000  X,  the  horizontal  sweep  of  the  beam  on  the 
record  screen  had  to  be  about  800  cm/sec.  This  gave  a  resolution  limit 
of  .04  cm  on  the  record  screen,  or  40  nm  on  the  sample.  (All  exposed 
areas  were  thus  40  nm  larger  in  the  horizontal  dimension  than  they 
would  have  been  without  this  resolution  problem.)  This  was  the  resolu¬ 
tion  limit  of  the  flying  spot  scanner  system.  The  use  of  a  record  scope 
tube  with  a  faster  phosphor  would  have  increased  the  resolution 
dramatically. 

This  resolution  of  the  flying  spot  scanner  was  not  the  factor  which 
limited  the  minimum  size  of  the  microbridges  produced.  As  mentioned 
above,  the  scattering  of  the  electrons  in  the  resist  limited  the 
minimum  size  of  the  exposed  pattern  to  about  one-fourth  of  the  resist 
thickness.^  In  all  samples  produced  this  scattering  limitation  was  larger 
than  the  flying  spot  scanner  limitation. 

The  flying  spot  scanner  system  worked  quite  well  in  practice;  Figure 
2.7  shows  a  structure  produced  using  this  system. 

2.4  Thin  Film  Deposition  and  Etching 

2.4.1  Thermal  evaporation 

Thin  films  of  Pb,  Sn,  In,  A1 ,  Cu,  Cr,  and  Ag  were  produced  in  a 
standard  evaporator.  The  source  to  substrate  distance  was  20  cm,  and  the 
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Figure 
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2.7  Uni  form- thickness  tin  microbridge  produced  with  EBL. 
The  tin  is  about  100  nm  thick. 
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thickness  was  monitored  with  an  Edwards  film  thickness  monitor.  The 
substrates  were  cooled  with  liquid  nitrogen  to  about  100°K.  Base 
pressures  before  evaporation  were  about  1  x  10‘6  Torr. 

A1  films  were  deposited  from  resistively  heated  tungsten  baskets, 
at  a  rate  of  about  1  nm/sec.  Cr  films  were  deposited  from  heated  Cr- 
plated  tungsten  rods,  at  a  rate  of  about  1  nm/sec.  Pb,  Sn,  In,  Au,  and 
Ag  were  deposited  from  heated  tantalum  boats,  at  rates  ranging  from  50 
nm/sec  for  lead  down  to  10  nm/sec  for  Cu. 

It  was  essential  to  cool  the  substrates  and  to  evaporate  the  films 
as  fast  as  possible  for  two  reasons:  it  was  necessary  have  grain  sizes 
less  than  100  nm,  and  the  resist  stencils  could  not  tolerate  much 
heating. 

All  of  the  films  produced  by  thermal  evaporation  were  patterned 
using  the  lift-off  method. 

2.4.2  Ion-beam  cleaning 

The  production  of  the  Pb-Cu-Pb  microbridges  required  that  the  Cu 
bridge  films  be  cleaned  immediately  before  the  Pb-bank  films  were 
deposited  in  order  to  remove  the  oxide  on  the  copper  and  the  resist 
residue  remaining  from  the  copper-stripe  fabrication  step.  This  cleaning 
was  accomplished  using  the  ion  mill  depicted  in  Figure  2.8.  The  ion 
source  was  purchased  from  Ion-Tech  Inc.;  the  electronics  were  constructed 
in-house  from  a  circuit  which  was  nearly  the  same  as  the  circuit  provided 
by  Ion-Tech.  The  ion  beam  had  a  diameter  of  2.5  cm  at  the  source,  and  a 
diameter  of  about  10  cm  at  the  substrates.  Thermal  evaporation  sources 
were  placed  alongside  the  ion  mill  so  that  films  could  be  deposited 
just  after  the  ion  cleaning  was  done  without  breaking  the  the  vacuum. 


The  Cu  films  were  ion-cleaned  for  150  sec  at  an  accelerating  voltage 

2 

of  600  volts  and  an  ion-current  density  of  about  0.2  ma/cm  at  the  sub¬ 
strate.  It  was  estimated  from  Cu  sputtering-yield  data11  that  about  10  nm 
of  Cu  was  removed  by  this  cleaning.  The  pressure  before  cleaning  was  less 
than  1  x  10“®  Torr;  the  argon  pressure  in  the  bell  jar  was  about  5  x  10"4 
Torr.  One  advantage  of  this  ion-cleaning  was  that  the  Pb  films  deposited 
after  ion-cleaning  were  much  more  adherent  than  Pb  films  deposited  without 
this  cleaning  step. 

Immediately  after  the  cleaning,  the  ion  mill  was  shut  off,  the  argon 
flow  was  stopped,  and  the  Pb  was  evaporated.  This  switch-over  from  ion 
cleaning  to  film  evaporation  took  about  10  sec;  this  delay  did  not  seem  to 
cause  any  noticeable  problems  in  the  Pb-Cu  interface.  Production  of 
Pb-Al-Pb  microbridges  in  the  same  fashion  was  attempted,  but  a  good  Pb-Al 
interface  could  not  be  obtained.  Either  the  A1  oxide  was  not  removed 
completely,  or  the  A1  was  re-oxidizing  during  the  switch-over  between 
cleaning  and  evaporation. 

2.4.3  Sputter  deposition  and  etching  of  Nb  films 

The  very  high  melting  temperature  of  niobium  prevented  the  thermal 
evaporation  of  niobium  in  the  available  apparatus,  so  niobium  films  were 
produced  by  sputter  deposition.  A  Varian  rf-diode  sputtering  system  was 
used.  This  system  had  a  base  pressure  of  2  x  10~7  Torr,  and  the  turbo¬ 
pump  used  did  not  have  to  be  throttled  during  sputtering;  the  contam¬ 
inating  residual  gases  were  thus  kept  to  a  minimum.  The  sputtering  gas 
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was  argon  at  a  pressure  of  10  Torr.  The  films  were  bias-sputter 
deposited  at  a  cathode  potential  of  1600  volts,  a  substrate  potential 
of  40  volts,  and  an  input  power  of  500  watts.  The  deposition  rate  was 


determined  by  interferometric  methods  to  be  0.65  nm/sec;  since  this  rate 
was  very  repeatable,  films  of  a  desired  thickness  were  made  by  sputtering 


for  a  pre-determined  time  interval.  Although  the  substrate  table  was  water- 
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cooled,  previous  measurements  indicated  that  during  deposition  the  sub¬ 
strate  temperature  could  rise  as  high  as  250°C.  The  films  produced  in 
this  system  were  of  high  quality:  films  100  nm  thick  had  a  Tc  of  9.1  K 
and  a  resistance  ratio  of  4;  films  35  nm  thick  had  a  T  of  5.6  K  and  a 
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resistance  ratio  of  3. 

A  few  attempts  were  made  to  use  the  lift-off  method  with  these 

films,  but  the  sputtering  process  proved  to  damage  the  resist  so  badly 

that  the  films  would  not  lift  off.  Consequently,  sputter  etching  was 

used  to  pattern  the  niobium  films  after  they  were  deposited.  Both  the 

el ectron-beam  resist  and  aluminum  were  used  as  etch  masks.  The  sputter- 
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etching  was  done  in  argon  at  5  x  10  Torr,  with  a  substrate  potential  of 
500  volts  and  an  input  power  of  100  watts.  The  etch  rate  was  about 
.086  nm/sec;  this  etch  rate  was  not  as  repoducible  as  the  deposition  rate, 
so  consecutive  etches,  each  etch  much  shorter  than  the  previous  etch, 
were  used  to  obtain  the  desired  removed  thickness  of  niobium. 

The  sputter  etching  could  also  be  used  to  clean  the  substrate 
and/or  previously  deposited  films  since  it  was  not  necessary  to  break  the 
vacuum  between  sputter  etching  and  sputter  deposition. 

2.5  Fabrication  Details  of  Specific  Samples 
2.5.1  Dayem  bridges 

The  first  samples  produced  in  this  research  effort  were  uniform 
thickness  (Dayem)  tin  bridges,  an  example  of  which  is  shown  in  Figure  2.7. 
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These  samples  were  produced  using  the  flying  spot  scanner  and  a  single 
tin  deposition.  Although  these  samples  were  quite  small  (500  nm  wide 
by  500  nm  long)  their  performance  was  only  mediocre.  Because  of  this 
poor  performance  they  were  studied  only  briefly. 

2.5.2  Indium  variable-thickness  bridges 

Indium  variable- thickness  microbridges  (Figure  2.9)  were  produced  in 
a  two-step  process,  both  steps  using  EBL  and  the  lift-off  method.  First, 
the  indium  bridge  line  was  laid  down  using  the  line-exposure  EBL  tech¬ 
nique;  these  bridge  lines  were  about  80  nm  thick.  Next,  the  indium  banks 
were  laid  down  using  the  flying  spot  scanner  technique;  these  films  were 
about  300  nm  thick.  These  microbridges  had  very  poor  performance,  probably 
due  to  very  bad  interfaces  between  the  bridge  films  and  the  banks  films; 
no  cleaning  steps  were  used  between  the  two  indium  depositions. 

2.5.3  Pb-Cu-Pb  microbridges 

The  most  successful  microbridges  produced  were  those  in  which  the 
birdge  was  made  from  copper  and  the  banks  were  made  from  lead  (Figure 
2.10).  The  fabrication  process  for  these  bridges  was  the  following: 

First,  two  registration  pads  were  put  on  the  substrate;  these  pads 
were  made  using  photolithography  and  the  lift-off  technique.  The  pads 
were  about  10  microns  wide  and  150  microns  apart;  they  consisted  of  80 
nm  of  silver  on  top  of  20  nm  of  chromium. 

Next,  the  copper  bridge  line  (200  microns  long,  about  0.2  micron 
wide)  was  placed  between  the  registration  pads  using  line-exposure  EEL 
and  the  lift-off  method.  The  copper  line  was  about  60  nm  thick. 

After  the  bridge  line  deposition,  the  resist  stencil  for  the  banks 
was  put  down.  This  pattern,  shown  in  Figure  2.3,  was  done  in  a  series  of 
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Figure  2.8  Diagram  of  the  ion  mill  -  evaporator  system. 
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Figure  2.9 


Indium  variable-thickness  microbridge.  The  bridge  is 
about  80  run  thick  and  the  banks  are  about  250  nm  thick. 


‘T  ■. 


Figure  2.10  Lead-copper  SNS  microbridge.  Although  the  banks  seem  to 
be  shorted  together,  electrical  measurements  showed  that 
they  are  not  (Pb-Cu-7). 
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three  consecutive  exposures;  the  first  exposure  was  done  with  the  flying 
spot  scanner  at  a  magnification  of  10,000  X.  The  second  and  third 
exposures  were  done  in  the  rectangle-exposure  mode,  at  mangifi cations  of 
1000  X  and  100  X,  respectively .  After  development,  the  copper  was  ion 
cleaned,  the  lead  banks  (about  150  nm  thick)  were  evaporated,  and  then 
the  unwanted  lead  was  lifted  off  by  dissolving  the  resist. 

The  final  fabrication  step  was  the  deposition  of  the  gross  elec¬ 
trodes.  A  narrow  (''-O.S  mm)  strip  of  aluminum  foil  was  taped  over  the 
bridge,  leaving  both  ends  of  the  lead  banks  sticking  out.  The  exposed 
lead  was  then  ion  cleaned  and  about  one  micron  of  lead  was  deposited;  the 
protective  aluminum  foil  strip  was  then  removed  and  the  sample  was  ready 
for  testing. 

2.5.4  Nb-Cu-Nb  and  Nb-AI-Nb  microbridges 

Microbridges  which  had  a  bridge  made  of  aluminum  or  copper  and  banks 
made  of  niobium  were  produced  using  a  combination  of  the  lift-off  method 
and  the  sputter-etch  method  to  pattern  the  films.  First,  registration 
pads  and  the  bridge  line  were  laid  down  in  a  fashion  similar  to  the  pro¬ 
cedure  used  for  the  Pb-Cu-Pb  structures;  in  these  cases  the  copper  films 
were  40  nm  thick  and  the  aluminum  films  were  40  nm  thick. 

Next,  the  bridge  was  cleaned  by  sputter-etching  and  a  niobium  film 
(65  nm  thick)  was  deposited  through  a  tantalum  foil  mask;  this  mask 
patterned  the  niobium  film  in  the  gross  electrode  structure  (Figure  2.2). 
Note  that  this  niobium  film  completely  covered  the  bridge  line. 

After  the  niobium  deposition,  the  rectangle-exposure  and  the  line- 
exposure  EBL  modes  were  used  to  produce  the  PMMA  etch  mask  on  top  of  the 
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niobium  film;  the  niobium  was  then  sputter  etched  to  form  the  bank 
geometry.  Care  had  to  be  taken  during  this  etch  in  order  that  all  the 
niobium  was  etched  out  of  the  gap  but  the  bridge  was  not  etched  through. 
Because  the  PMMA  sputter  etched  three  times  faster  than  the  niobium, 
and  because  the  resist  could  only  be  about  300  nm  thick,  the  niobium  banks 
had  to  be  less  than  100  nm  thick.  In  addition,  the  sputter-etching 
process  hardened  the  resist  making  it  difficult  to  remove;  fortunately 
this  unsightly  resist  residue  did  not  affect  the  microbridge's  per¬ 
formance  . 

Finally,  thick  1  micron)  gross  niobium  electrodes  were  deposited 
using  a  technique  similar  to  that  used  for  the  final  gross  electrode 
fabrication  of  the  Pb-Cu-Pb  microbridges,  except  that  sputter  etching 
and  sputter  deposition  were  used.  An  example  of  the  final  result  is 
shown  in  Figure  2.11 . 

The  performance  of  these  Nb-Cu-Nb  and  Nb-Al-Nb  microbridges  will 
be  discussed  in  Chapter  5.  In  general,  the  bridges  had  long  effective 
lengths  as  a  result  of  short  electron  mean  free  paths  in  the  bridges. 

These  samples  had  relatively  small  ICR  products,  and  the  Nb-Al-Nb 
samples  had  poor  response  to  microwave  radiation. 

2.5.5  Nb-Nb-Nb  microbridges 

Variable-thickness  microbridges  in  which  both  the  bridge  and  the 
banks  were  niobium  were  produced  by  depositing  a  niobium  film  and  then 
using  sputter-etching  to  produce  both  the  bridge  and  the  banks.  First, 
a  niobium  film  (10C  nm  thick)  was  deposited  thorugh  a  tantalum  mask. 

Next,  a  thin  aluminum  line  (40  nm  thick)  was  laid  over  the  top  using 
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Figure  2.11  Niobium-copper  SNS  microbridge. 
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line-exposure  EBL  and  lift-off.  This  line  was  used  as  an  etch  mask  to 
form  a  raised  line  in  the  niobium  by  sputter-etching  away  about  35  nm 
of  the  niobium.  The  aluminum  was  then  removed  in  NaOH. 

Next,  the  gap  for  the  banks  was  sputter  etched  in  a  fashion  similar 
to  that  used  to  produce  the  gap  in  the  niobium  in  the  Nb-Cu-Nb  structures. 
Again  care  had  to  be  taken  to  ensure  that  the  etch  was  long  enough  to 
remove  all  the  niobium  in  the  gap,  but  not  so  long  that  the  bridge 
(corresponding  to  the  raised  portion  of  the  niobium)  was  not  etched 
completely  away.  Since  both  aluminum  and  PW1A  etch  much  faster  than 
niobium,  the  original  niobium  films  had  to  be  less  than  100  nm  thick.  This 
meant  that  the  bridge  and  banks  were  relatively  thin  compared  to  the 
length  and  width  of  the  bridge,  resulting  in  poor  performance  for  these 
microbridges.  Figure  2.12  shows  one  of  the  bridges  of  this  type,  and 
Figure  2.13  shows  a  photograph  of  one  of  these  samples. 


Figure  2.12  Niobium  variable- thickness  microbridge.  The  bridge  Is 
about  25  nm  thick  and  the  banks  are  about  60  nm  thick 
(Nb37-2). 


Figure  2.13  Photograph  of  the  sample  shown  ir  Figure  2.12. 


Chapter  3 

MEASUREMENT  APPARATUS  AND  TECHNIQUES 

The  measurement  of  the  electrical  properties  of  the  microbridges 
required  an  apparatus  to  cool  the  samples  and  electronics  with  which  to 
make  the  actual  measurements.  Since  our  laboratory  has  had  many  years 
of  experience  in  cryogenic  technology,  obtaining  a  cryostat  in  which  to 
cool  the  samples  was  relatively  straightforward.  On  the  other  hand,  the 
electronics  used  in  this  work  was  somewhat  different  from  that  used 
previously  in  our  laboratory,  primarily  due  to  the  fact  that  the  quality 
of  the  electronic  equipment  available  was  constantly  improving.  In  light 
of  these  facts,  the  cryogenic  apparatus  will  be  dealt  with  in  a  general 
manner,  while  the  electronics  will  be  discussed  in  more  detail. 

3.1  The  Cryostat 

The  cryostat  used  in  this  work  was  originally  built  by  Dr.  James 
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Lukens,  but  many  changes  have  been  made  to  it  since  his  work.  Jackel 
has  described  the  construction  and  operation  of  this  cryostat  in  detail. 
Since  we  made  no  changes  in  the  thermal  construction  and  operation  of 
this  apparatus,  the  reader  is  referred  to  Jackel  for  a  description  of 
these  aspects.  Almost  all  of  our  alterations  concerned  the  electrical 
wiring  of  the  cryostat,  and  these  aspects  will  be  considered  here. 

Figure  3.1  shows  a  schematic  diagram  of  the  lower  part  of  the 
cryostat;  only  the  electrical  features  are  shown  in  detail.  There  were 
two  distinct  sets  of  electrical  connections  to  be  made;  the  first  was 
the  dc  current  and  dc  voltage  connection,  and  the  second  was  the 
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Figure  3. 


1  Schematic  diagram  of  the  lower  portion  of  the  cryostat. 
For  clarity  most  of  the  temperature  control  components 
have  not  been  included. 
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Figure  3.1 


microwave  connection.  In  implementing  both  of  these  sets  of  connections 
it  was  necessary  to  preserve  the  vacuum  integrity  of  the  vacuum  can  and 
the  exchange-gas  can,  as  well  as  preserving  the  weakness  of  the  thermal 
connection  between  these  two  cans. 

3.1.1  Sample  holder 

Figure  3.2  shows  a  schematic  diagram  of  the  sample  holder.  The 
sample  was  clamped  between  the  phenolic  block  and  a  piece  of  transparent 
epoxy,  with  the  electrical  connections  made  by  spring-loaded  indium- 
coated  brass  pins.  The  epoxy  piece  was  made  transparent  so  that  the 
sample  could  be  observed  during  mounting  to  ensure  that  the  pins  made 
contact  at  the  correction  locations.  This  sample  holder  worked  very  well; 
there  was  seldom  a  loss  of  electrical  contact  during  thermal  cycling,  and 
the  sample  could  be  removed  and  then  remounted  at  a  ~ater  date  for  addi¬ 
tional  testing.  A  photograph  of  this  holder  is  shown  in  Figure  3.3. 

3.1.2  DC  I -V  wiring 

The  dc  I-V  wiring  can  be  seen  in  Figure  3.1.  Because  the  micro¬ 
bridges  were  quite  small,  they  had  to  be  protected  against  burn-out 
caused  by  high-frequency  induced  currents.  To  protect  the  sample,  a  low- 
pass  filter  was  placed  in  series  with  the  I-V  leads.  This  filter  is  shown 
in  Figure  3.4.  These  filters  effectively  shorted  out  the  sample  at  fre¬ 
quencies  above  about  1  kHz,  while  allowing  dc  measurements  to  be  per¬ 
formed  below  about  1  kHz;  this  upper  limit  on  the  measurement  frequency 
was  not  a  hindrance  to  the  measurements. 


In  order  to  protect  the  sample  from  dc  current  surges  when  the 
electronics  was  connected,  a  room- temperature  shorting  box  was  put  at  the 
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Figure  3.2  Schematic  diagram  of  the  sample  holder. 
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Figure  3.3  Photograph  of  the  sample  holder. 
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top  of  the  cryostat.  This  box  had  a  switch  which  could  short  the  I  leads 
together  and  the  V  leads  together  while  the  electronics  was  being 
connected,  and  then  isolate  the  leads  while  measurements  were  made.  This 
box  is  depicted  in  Figure  3.4. 

To  obtain  the  required  thermal  isolation,  the  I-V  leads  were  brass 
from  room- temperature  down  to  the  4.2K  filters,  niobium  from  the  filters 
down  to  the  shielding  can,  and  copper  from  there  to  the  sample  holder. 

The  wires  were  fed  through  the  vacuum  can  and  the  exchange-gas  can  via 
vacuum-tight  epoxy  feed-throughs. 

3.1.3  Microwave  connection 

The  sample  could  be  fed  by  a  microwave  current  by  inductively 
coupling  the  sample  to  a  microwave  coaxial  line  via  a  loop  which  terminated 
the  line.  Although  no  serious  attempt  was  made  to  optimize  this  coupling, 
sufficient  coupling  was  produced  for  our  purposes.  The  coaxial  line 
had  a  silver-plated  steel  inner  conductor  and  a  stainless-steel  outer 
conductor  from  room  temperature  down  to  the  shielding  can.  Although  this 
coax  had  a  fairly  high  attenuation,  it  was  necessary  in  o^der  to  achieve 
the  desi”ed  thermal  isolation.  From  the  shielding  can  to  the  sample 
both  the  inner  and  outer  conductors  were  made  of  copper. 

3.1.4  Measurements  below  1.5K 

This  cryostat  was  useful  down  to  about  1.5K.  In  order  to  make  dc 
I-V  measurements  down  to  about  0.5K,  a  cryostat  which  was  cooled  by 
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pumped  He  was  used;  this  cryostat  was  cheerfully  provided  by  Alan 


Kleinsasser. 


3.2  Electronics 


Several  different  configurations  of  electronics  were  used  since 
different  types  of  measurement  were  made.  In  general,  since  the  voltages 
and  currents  of  interest  were  small  (as  low  as  100  nA  and  100  nV),  it 
was  necessary  to  use  the  lowest-noise  electronics  available.  But,  since 
the  available  electronics  quite  often  produced  a  small  signal-to-noise 
ratio,  it  was  often  necessary  to  average  the  measurements  for  long 
periods  of  time. 

Extreme  care  was  taken  to  minimize  the  effects  of  an  electrically 
noisy  environment.  The  dewar  was  in  a  shielded  enclosure,  and  most  of  the 
electronics  were  placed,  along  with  the  dewar,  in  a  shielded  room. 

Digital  electronics  had  to  be  kept  outside  the  shielded  room  due  to  their 
excessive  noise  radiation.  In  order  to  reduce  60  Hz  noise,  which  the 
shielded  room  could  not  screen  out,  all  wires  were  either  twisted  pairs 
or  triaxial  shielded  cable.  Ground  loops  were  very  carefully  eliminated 
by  having  all  circuits  grounded  at  only  one  point. 

3.2.1  Basic  I-V  measurement  system 

Figure  3.5  shows  the  basic  configuration  for  dc  I-V  measurements  or 
dV/dl  measurements.  The  sample  was  current  biased  by  a  current  source 
(shown  in  Figure  3.6).  This  current  source  produced  a  current  ramp,  or 
it  converted  a  voltage  ramp  from  a  triangular  oscillator  into  a  current 
ramp.  The  current  was  measured  by  a  current-monitoring  resistor  (CMR); 
either  a  resistor  in  the  current  source  or  a  resistor  in  the  shorting 
box  was  used  for  a  CMR.  The  voltage  across  the  CMR  was  fed  into  the  input 
of  a  PAR  113  low-noise  differential  amplifier  and  the  output  of  the  113 


Figure  3.5  Schematic  diagram  of  the  basic  measurement  electronics. 
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Figure  3.6  Schematic  diagram  of  the  current  source. 


went  into  the  horizontal  axis  of  an  x-y  recorder. 


The  voltage  across  the  sample  was  fed  directly  into  a  PAR  113  pre¬ 
amplifier,  and  the  output  of  this  113  went  into  the  vertical  axis  of  the 
x-y  recorder. 

For  dV/dl  measurements,  a  small  ac  voltage  (from  the  reference  of 
a  PAR  124  lock-in  amplifier)  was  also  fed  into  the  current  source, 
resulting  in  a  small  ac  component  to  the  sample  bias  current.  (Usually, 
an  ac  frequency  of  1  kHz  was  used,  with  ar,  ac  magnitude  of  about  one 
percent  of  the  peak  bias  current.)  The  output  of  the  voltage-measuring 
PAR  113  was  now  fed  into  the  vertical  axis  of  the  x-y  recorder.  The 
signal  from  the  voltage  113  and  the  dV/dl  were  plotted  on  the  same  piece 
of  paper,  resulting  in  curves  like  the  ones  shown  in  Figure  3.7. 

The  basic  problem  with  this  configuration  was  that  if  the  signal-to- 
noise  level  was  small,  very  slow  ramps  were  necessary.  Not  only  was  it 
difficult  to  generate  very  slow  ramps,  but  the  dc  drift  in  the  pre¬ 
amplifiers  put  a  lower  limit  on  the  speed  of  the  ramp. 

3.2.2  Digital  data  acquisition 

Because  of  the  problem  with  dc  drifts  in  the  preamplifiers,  a 
digital  data  acquisition  system  was  often  used.  This  system  is  shown 
in  Figure  3.8.  In  this  setup,  the  dc  current-ramp  biasing  the  sample 
was  repeated  many  times,  the  corresponding  voltage  or  dV/dl  curves  for 
each  separate  ramp  were  added  together  digitally  with  a  Northern  signal 
averager.  If  the  time  for  one  complete  ramp  was  less  than  the 
characteristic  dc  drift  time  of  the  preamp,  the  dc  drift  of  tPie  preamp 
would  appear  as  a  constant  offset  in  the  averaged  curve.  Thus  the  system 
allowed  the  data  to  be  averaged  over  long  periods  of  time  without  the 
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Figure  3.7  Typical  high-voltage  I-V  and  dV/dl  curves  generated  with 
the  basic  electronics  configuration.  The  curves  are 
symmetric  about  the  I  =  0  point.  Both  curves  have  the 
same  current  (horizontal)  scale;  the  vertical  scale  for 
dV/dl  is  in  arbitrary  units. 
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dc  drift  changing  the  shape  of  the  curves. 

In  this  configuration  the  sample  bias  current  was  not  measured 
directly,  but  was  measured  by  assuming  the  current  to  be  proportional 
to  the  amount  of  time  since  the  sync  pulse  of  the  oscillator  producing 
the  ramp. 

For  most  of  the  data  taken  in  this  work  the  digitally-averaged  data 
were  far  superior  to  the  data  produced  using  the  more  basic  system;  in 
fact,  nearly  all  of  the  experimental  data  presented  in  this  thesis  was 
taken  using  this  system.  This  digital  configuration  was  not  without  its 
problems,  however.  The  most  serious  disadvantage  was  that  the  zero  of 
the  vertical  axis  was  not  known;  in  fact,  the  offset  from  zero  varied 
from  one  curve  to  the  next.  This  prevented  using  this  system  to  deter¬ 
mine  the  temperature  dependence  of  the  magnitude  of  dV/dl;  the  more 
basic  system  had  to  be  used  in  this  case. 

Another  serious  disadvantage  was  that  hysteresis  in  the  data  was  not 
easily  observed.  To  observe  hysteresis,  it  was  necessary  to  separately 
average  the  curves  for  increasing  bias  current  and  the  curves  for 
decreasing  bias  current,  and  then  compare  the  averaged  curves.  This  was 
something  of  a  nuisance. 

The  digital  system  was  not  as  accurate  as  the  more  basic  system.  The 
non-linearity  in  the  A-D  converter  of  the  signal  averager  was  about  one 
percent.  In  addition,  the  gain  of  the  signal  averager  was  not  accurately 
known;  the  peak-to-peak  voltage  of  the  output  of  the  preamp  had  to  be 
measured  on  an  oscilloscope  to  determine  the  gain.  A  similar  measurement 
had  to  be  done  to  measure  the  current  gain. 
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3.2.3  Digital  data  processing 

One  advantage  of  the  digital  data  acquisition  system  was  that  the  data 
could  not  only  be  plotted  on  an  x-y  recorder,  but  it  could  also  be  sent 
directly  to  a  Prime  computer.  This  allowed  the  data  to  be  manipulated 
by  the  computer  before  it  was  plotted  out.  For  example,  in  measuring  the 
dV/dl  data  the  direct  result  was  a  curve  of  dV/dl  versus  current.  But 
since  the  signal  averager  could  simultaneously  acquire  both  the  dV/dl 
versus  current  curve  and  the  voltage  versus  current  curve,  it  was  easy  to 
program  the  computer  to  plot  dV/dl  versus  voltage.  The  dV/dl  versus 
voltage  curves  presented  in  Chapter  4  were  obtained  in  this  way. 

It  was  also  possible  to  derive  dV/dl  versus  voltage  curves  from 
voltage  versus  current  curves  with  the  computer.  An  example  of  this  is 
shown  in  Figure  3.9.  In  this  way  it  was  possible  to  make  precise  measure¬ 
ments  of  the  voltages  Vinf^-j ecti0n  and  Vlower  discussed  in  Chapter  4. 

Having  the  data  stored  in  the  computer  also  made  it  possible  to 
replot  the  data  at  any  time,  at  any  scale;  if  necessary,  very  small 
regions  of  a  curve  could  be  blown  up,  or  two  different  curves  could  be 
compared.  This  was  very  useful  in  analyzing  the  data,  since  it  was 
possiole  to  look  at  features  of  the  data  which  were  not  necessarily 
considered  noteworthy  when  the  data  was  taken. 

3.2.4  Microwave  electronics 

Figure  3.10  shows  a  schematic  of  the  electronics  used  to  bias  the 
sample  with  an  rf  current.  For  frequencies  less  than  9  GHz,  a  TWT 
oscillator  was  used  as  a  source  of  microwaves.  In  this  case  the  post- 
directional  -coupler  attenuator  was  held  constant,  and  the  pre-directional- 
coupler  attenuator  was  changed  to  vary  the  magnitude  of  the  applied  rf. 


The  power  (or,  equivalently,  the  magnitude  of  the  rf  current)  was  deter¬ 
mined  by  observing  the  power  meter.  For  these  frequencies  all  of  the 
components  were  connected  by  coaxial  lines. 

For  frequencies  above  9  GHz  a  klystron  was  used  as  a  microwave 
source.  In  this  case  the  magnitude  of  the  applied  rf  was  varied  by 
changing  the  post-direction-coupler  attenuator,  and  the  power  applied  to 
the  sample  was  measured  by  noting  the  attenuator  setting.  The  pre- 
directional-coupler  attenuator  was  used  to  hold  the  power  to  the  post- 
directional  -coupler  attenuator  constant,  as  measured  by  the  power  meter. 
For  these  frequencies  all  of  the  components  up  to  the  final  attenuator 
were  connected  by  rectangular  waveguide,  in  order  to  minimize  excess 
attentation.  A  waveguide-to-coax  adaptor  was  inserted  between  the  final 
attenuator  and  the  cryostat's  coaxial  line. 

For  all  frequencies,  it  was  necessary  to  tune  the  frequency  for 
maximum  rf  power  coupled  into  the  sample  because  of  the  strong  mismatches 
present  in  the  cryostat's  coaxial  lines  and  the  mismatch  between  the 
terminating  loop  and  the  sample. 


Chapter  4 

RESULTS  ON  Pb-Cu-Pb  JUNCTIONS 

The  measurements  of  the  electrical  properties  of  the  SNS  junctions 
produced  in  this  work  were  multi-faceted;  in  light  of  this,  the  results 
are  arranged  here  in  separate  sections,  with  the  appropriate  theoretical 
and  experimental  backgrounds  presented  with  each  separate  aspect  of  the 
observed  results.  Tnis  division  of  results  is  somewhat  artificial, 
however,  in  that  the  behavior  seen  in  one  type  of  measurement  was  often 
very  useful  in  understanding  the  results  of  other  measurements.  Some 
of  the  results  were  fairly  simple  in  nature,  and  the  understanding  of 
these  results  was  relatively  simple.  Other  results  were  complex,  and 
the  understanding  of  these  results  was  correspondingly  complex.  A  few 
of  the  measurements  showed  quite  unexpected  effects,  and  some  of  these 
effects  are  still  poorly  understood. 

4.1  Temperature  Dependence  of  the  Critical  Current 

4.1.1  Theory 

The  existence  of  supercurrents  in  an  SNS  microbridge  at  first 

thought  seems  surprising  since  the  bridge  is  made  from  a  normal  metal. 

After  all,  the  electron-phonon  interaction  responsible  for  super- 

13 

conductivity  has  a  range  of  only  a  few  interatomic  distances  and 
allowing  for  retardation  effects  only  increases  the  interaction  length 
to  about  20  nm;  the  microbridges  produced  in  this  work  had  normal - 
metal  bridges  200  nm  long  or  longer.  The  answer,  of  course,  is  that 
there  is  still  a  coupling  between  the  superconducting  banks.  This 
coupling  is  due  to  the  fact  that  a  superconducting  bank  imposes  on  the 
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normal  bridge  a  change  in  the  density  of  states.  It  is  important  to  note 

that  this  change  in  the  density  of  states  does  not  include  an  induced 

14 

energy  gap  in  the  normal  metal. 

The  bank  also  induces  a  correlation  in  the  electron  distribution 
in  the  normal  metal.  In  the  superconductor  some  of  the  electrons  are 
paired  in  states  of  opposite  momentum;  this  pairing  between  electrons 
extends  into  the  normal  metal  (proximity  effect).  This  correlation  decays 
exponentially  with  distance  into  the  normal  metal,  with  the  decay  length 
given  by 


5(T) 


M 

6irkT 


1/2 


(vf  =  fermi  velocity  in  the  normal  metal,  A  =  electron  mean  free  path  in 
the  bridge).  In  spite  of  this  decay,  the  opposite  superconductor  still 
"sees"  some  of  this  correlation;  i.e.,  the  two  superconductors  are 
coupled  together.  It  is  this  coupling  which  gives  rise  to  supercurrents 
in  the  normal  metal.  Figure  4.1  illustrates  this  idea. 

Putting  this  electron-correlation  picture  of  SNS  structures  on  a 
firm  theoretical  footing  has  been  the  subject  of  a  great  deal  of  work. 

The  calculation  of  the  density  of  states  near  the  interface  will  be  dis¬ 
cussed  later  in  this  chapter.  The  calculation  of  the  critical  current  has 
been  carried  out  using  several  different  methods. 

Microscopic  treatments  for  clean  systems  (mean  free  path  much  greater 

than  the  junction  length)  were  performed  by  Kulik^,  Ishii^,  Bardeen 
17  18 

and  Johnson  ,  and  Bezuglyi  et  aj..  Ishii  restricted  himself  to  the 
T*0  K  case.  The  remainder  of  these  authors  predicted  that  Ic  should  be 
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Figure  4.1  Schematic  diagram  of  the  spatial  variation  of  the 
paired-electron  density  ( | F | )  in  a  SNS  microbridge. 
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proportional  to  e"'". 

19 

Kulik  and  Mltsai  microscopically  calculated  I  versus  T  for  a 

moderately  dirty  (£  «  A  «  L)  junction  and  found  I  ^exp-[(|^  +  Jr)]. 

1  c  ^o  h 

DeGennes  used  a  Ginzburg-Landau  treatment  for  the  banks  and  a 

2 

microscopic  treatment  for  the  bridge  and  found  I  ^  (1  -  T/T  )  ,  where 

Tc  is  the  transition  temperature  of  the  banks. 

20 

Fink  ,  using  Ginzburg-Landau  theory  for  both  bridge  and  banks,  found 

Ic  n,  (1  -  T/Tc)n  where  1  <  n  <  2.  The  calculations  of  both  deGennes  and 

Fink  are  restricted  to  the  dirty  limit  near  T  . 

21 

Likharev  has  done  a  microscopic  calculation  in  the  dirty  limit 

22 

using  the  equations  of  Usadel .  Since  this  calculation  has  been  very 
useful  in  this  work,  it  will  be  explained  here  in  some  detail.  In  this 
theory  the  correlation  picture  of  superconductivity  is  put  on  a  rigorous 
footing  by  defining  a  correlation  parameter 


F(r)  =  <*+(r)*+(r)> 


where  ip(r)  is  the  one-electron  operator.  F(r)  is  thus  the  probability 
of  finding  two  electrons  in  the  condensed  (correlated)  state  at  the 
point  r.  Usadel  developed  a  set  of  diffusion-like  equations  for  deter¬ 
mining  F.  He  used  the  expression  for  F  obtained  in  the  temperature- 
dependent  Green's  function  method;  thus  he  derived  an  infinite  set  of 
equations.  For  a  one-dimensional  system  like  the  microbridges  produced 
here,  Usadel 's  equations  reduce  to 


2uF(»,x)  -  D^[l  -  lF(a.,x)  l2]^ 

+  1  - F^x) -  dlFiu>axl[5\<2Mx)[l  -  |F(u.,x)12]>t 

2  [1  -  |F(u,x)!ZP  dx  ) 


■uuMHitaMiaAiii 
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v 

where  D  s  -j-  is  the  diffusion  constant  (v^  *  fermi  velocity,  A  *  electron 
mean  free  path),  u>  =  Matsubara  frequency  =  (2n+l)TrkT  (n  a  non-negative 
integer),  and 

A(x)  =  XN(0)27rkT  J  <F(u>,x)> 
oj>0 

(N(0)  =  density  of  states  at  fermi  surface,  X  =  interaction  constant). 

Note  that  all  of  these  equations  are  coupled  together  via  the  expression 
for  A;  this  coupling  makes  the  general  solution  to  these  equations 

s 

extremely  difficult  to  obtain.  But  for  a  normal  metal  X  =  0,  a  fact 
which  uncouples  the  equations  and  greatly  simplifies  finding  the  solution. 
Usadel  also  found  the  expression  for  the  supercurrent  density  to  be 

js(x)  =  2ieN(0KTD  F*(w,x)  ^  F(u,x)  -  F(u>,x)[^  F(w,x)]* 

It  is  important  to  note  that  these  equations  of  Usadel  are  valid  only  in 
the  dirty  limit. 

Following  Likharev  we  obtained  the  solution  to  these  equations  by 
numerical  integration;  these  calculations  provided  an  ICR  versus  T  relation 
which  depends  on  the  length  of  the  bridge.  Figure  4.2  shows  the  results 
of  Likharev's  calculations.  We  duplicated  Likharev's  calculations  in 
order  to  get  precise  theoretical  curves  with  which  to  fit  our  data;  the 
details  of  these  numerical  calculations  are  given  in  the  Appendix.  In 
these  calculations  it  was  necessary  to  impose  boundary  conditions  on  F 
at  the  superconductor-normal  metal  interface.  We  assumed  that  |F|  at 
the  boundary  had  the  same  value  as  it  did  in  the  bulk  of  the  super¬ 
conductor;  i.e. 
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Ip  |  s  ..  a.LiJ 
“  banks  [<d2  +  A(T)^]Js 

where  A(T)  is  the  temperature-dependent  gap  in  the  bulk  of  the  super¬ 
conductor.  This  assumption  is  not  valid  near  T  where  the  gap  in  the 
superconductor  is  depressed  near  the  bridge  due  to  the  proximity  of  the 
normal-metal  bridge.  DeGennes7  has  calculated  this  depression,  but 
application  of  his  boundary  condition  to  this  problem  was  not  attempted 
since  it  requires  a  better  knowledge  of  the  nature  of  the  bank-bridge 
interface  than  we  had  available. 

These  equations  for  F  are  somewhat  simplified  if  lengths  are  scaled 
in  units  of  the  decay  length  (coherence  length) 


Thus  a  bridge  of  physical  length  L  can  be  said  to  have  a  reduced  length 

i  =  yTt — - - r  ;  i.e.»  the  length  of  the  bridge  in  units  of  the  value 

^'e-banks' 

of  the  decay  length  at  the  transition  temperature  of  the  banks. 

4.1.2  Previous  experimental  work 

There  have  been  two  previous  measurements  made  of  the  versus  T 

23 

relation  for  SNS  junctions.  Clarke  made  measurements  on  layered  (as 

opposed  to  planar)  Pb-Cu-Pb  junctions.  His  junctions  had  a  cross- 

sectional  area  of  about  .04  mm  .  He  found  I  (1  -  T/Tc)  between  about 

6.5  K  and  7.0  K,  in  agreement  with  the  theory  of  deGennes.7  Below  this 

temperature  range  the  critical  current  was  limited  by  its  own  magnetic 

field  due  to  the  relatively  large  sire  of  the  junctions. 

24 

Shepherd  also  measured  layered  Pb-Cu-Pb  junctions;  his  samples 

2  -T 

had  an  area  of  about  .1  mm  .  He  found  that  I  n.  e  for  3K  £  T  <_  5K. 


iittk 


Since  his  samples  were  in  the  clean  limit,  this  result  is  inagreement  with 
the  microscopic  calculations  for  clean  junctions.  Below  5K  the  critical 
current  of  his  samples  was  self-limited. 


The  most  significant  difference  between  these  previous  experiments 
and  the  measurements  of  this  work  is  in  the  size  of  the  junctions. 
DeGennes^  has  calculated  that  for  a  400  nm  long  junction  the  penetration 
depth  is  on  the  order  of  1  mm.  Since  our  samples  had  widths  and  thick¬ 
nesses  less  than  .001  mm,  the  critical  current  was  uniform  across  the 
entire  cross  section  of  the  bridge  at  all  temperatures;  i.e.,  there  was 
no  self-limiting  of  the  critical  current,  and  meaningful  comparisons  with 
the  theory  were  made  down  to  the  lowest  temperature  achievable  (0.5K). 

4.1.3  Experimental  results 

Figure  4.3  shows  experimental  critical  current  versus  reduced 
temperature  (Tc  «  7.26K)  data  for  several  samples.  (The  critical  currents 
for  each  sample  have  been  multiplied  by  the  normal -state  resistance  RQ 
for  that  sample;  the  determination  of  RQ  will  be  considered  in  the 
low-voltage  I-V  section  of  this  chapter.)  The  current  is  plotted  on  a 
logarithmic  scale,  and  it  is  significant  to  note  that  the  data  for  the 
shorter  samples  cannot  be  fit  by  a  straight  line  at  low  temperatures. 

This  lack  of  a  straight-line  fit  implies  that  the  critical  current  at  low 
temperatures  is  not  proportional  to  a  power  of  T.  (This  is  not  in  con¬ 
flict  with  those  theories  which  predict  an  exponential  behavior  since 
those  theories  were  done  in  the  clean  limit,  while  our  samples  are 
clearly  in  the  dirty  limit.) 

21 

Figure  4.4  shows  a  fit  of  our  Likharev-like  calculations  to  the 
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Figure  4. 


3  Experimental  data  on  the  temperature  dependence  of  I  R 
for  four  Pb-Cu-Pb  microbridges  (Tc  =  7.26  K). 
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data  for  Pb-Cu-13.  Because  the  mean  free  path  A  and  the  physical  length 
L  were  not  accurately  known  for  our  samples,  the  reduced  length  l  could 
not  be  directly  calculated.  Instead,  i  was  determined  by  using  it  as  a 
fitting  parameter  in  fitting  the  experimentally  determined  I£R  versus  T 
curves  with  the  theoretical  curves.  We  found  that  for  the  shorter 
samples  it  was  necessary  to  increase  the  resistance  slightly  in  order  to 
fit  the  shape  of  the  theoretical  curves  to  the  experimental  data.  This 
can  be  seen  in  Figure  4.5.  The  calculation  for  l  =  8.9  fits  the  magnitude 
of  the  data  best,  but  the  a  •  8.2  fit  matches  the  shape  best.  There  are 
two  possible  explanations  for  our  data  to  have  a  lower  I  R  than  predicted. 

v 

The  measured  resistance  might  be  too  low  since  the  portion  of  the  bridge 
very  close  to  the  banks  might  not  contribute  to  the  measured  resistance 
since  the  density  of  pairs  is  very  high  there;  this  effect  would  be  more 
pronounced  in  shorter  bridges. 

The  other  possible  explanation  is  that  since  the  bridge  and  the  bank 
have  different  densities  of  states,  the  assumption  that  |F|  is  continuous 
across  the  interface  might  be  slightly  wrong. 

We  have  chosen  to  adjust  the  resistance  so  that  the  data  fits  the 
theory  at  low  temperatures.  We  have  thus  taken  the  l  value  determined 
in  this  way  to  be  the  correct  reduced  length  for  the  bridge.  In  any 
case,  the  different  possible  l  values  only  differ  by  less  than  10%. 

Table  1  gives  the  reduced  lengths  and  the  corrected  resistances 
used  to  obtain  them  for  some  of  the  samples  measured. 

The  success  of  this  theoretical  calculation  in  fitting  the  experi¬ 
mental  data  has  two  important  consequences.  First,  it  provides  a  good 
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Figure  4.5  Fit  of  the  theoretical  ICR  versus  T  relation  for  a  bridge 
of  length  £.  =  8.9  to  the  data  of  Pb-Cu-13.  The  data  has 
not  been  rescaled. 
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test  of  the  validity  of  Usadel's  equations.  This  appears  to  be  the  first 
experimental  verification  of  these  equations.  Second,  it  gives  reduced 
lengths  £  for  the  measured  samples.  Knowledge  of  these  lengths  has 
proven  very  useful  in  understanding  the  low-voltage  behavior  of  these 
samples. 

It  can  be  seen  in  Figure  4.4  that  the  Likharev  calculation  does  not 

fit  the  data  at  high  temperatures.  This  is  to  be  expected  since  this 

theory  ignores  the  effect  of  the  bridge  in  reducing  the  order  parameter 

in  the  banks,  an  effect  which  is  non-negligible  close  to  T  Close  to 

Tc,  however,  the  theory  of  deGennes^  should  be  valid.  Figure  4.6  shows 

o 

a  plot  of  I  versus  (1  -  T/T  )  close  to  T  for  Pb-Cu-13.  The  data  can 

C  C  w 

be  fit  by  a  straight  line  in  agreement  with  the  theory  of  deGennes. 
Similar  plots  for  the  other  samples  are  not  possible  since  the  critical 
currents  for  those  samples  were  too  small  to  be  observed  close  to  Tc  of 
the  banks. 

4.1.4  Current-phase  relation 

We  have  also  measured  the  current-phase  relation  for  these  samples. 
That  is,  the  supercurrent  is  dependent  on  the  phase  difference  between 
the  order  parameters  of  the  opposing  superconducting  banks, 

Ic  *  IQf (A<j>) ,  and  we  have  measured  the  functional  form  of  f.  The  details 

3 

of  these  measurements  has  been  given  by  Jackel .  We  found  that  for  all 
accessible  temperatures  f ( A4>)  =  sinA$. 

All  of  the  microscopic  theories  discussed  above  have  predicted  that 
for  very  short  SNS  samples  the  current-phase  relation  should  be  non- 
sinusoidal  at  low  temperatures .  But  for  samples  as  long  as  the  ones  used 
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Figure  4.6  Critical  current  of  Pb-Cu-13  plotted  versus  (1-T/T  ) 
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for  T/T  >  0.9. 
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in  this  work,  a  sinusoidal  current-phase  relation  has  been  predicted. 

The  details  of  this  prediction  are  given  in  the  Appendix  in  the  discussion 
on  the  integration  of  Usadel's  equations. 

4.2  Low  Voltage  I-V  Characteristics 

4.2.1  Experimental  results 

Figures  4.7  through  4.9  show  the  low-voltage  critical  current  versus 
voltage  (I-V)  curves  for  some  of  the  samples  measured.  All  of  these  I-V 
characteristics  have  certain  features  in  common.  For  almost  all  tempera¬ 
tures  below  the  transition  temperature  of  the  banks  there  is  a  critical 
current;  above  this  critical  current  a  voltage  appears  across  the  sample. 
Between  two  characteristic  voltages  V1/NWQW  and  V1M,M„  the  I-V  curve  is  a 
straight  line.  (Vuprar  is  not  shown  in  these  figures  but  it  is  shown  in 
later  figures  and  its  magnitude  for  various  samples  is  found  in  Table  1.) 
For  a  given  sample  the  value  of  the  slope  of  this  straight  line  is  inde¬ 
pendent  of  temperature,  except  for  very  close  to  T  .  (Near  Tc  the  slope 
of  this  line  increases  with  increasing  temperature.  In  Part  4  of  this 
chapter  it  will  be  shown  that  this  increase  in  slope  with  temperature  is 
most  likely  due  to  an  extra  dissipation  occurring  ir.  the  banks.)  We 
have  chosen  the  low  temperature  value  of  this  slope  to  be  the  normal- 
state  resistance  RQ  of  the  bridge.  RQ  was  the  resistance  required  in 
Part  1  to  fit  the  temperature  dependence  of  the  critical  current  to 
theory . 

The  character  of  the  I-V  curves  above  the  critical  current  but 
below  VUpper  can  be  classified  by  three  temperature  regimes.  At  high 
temperatures  the  I-V  curve  is  continuous,  rising  smoothly  from  the 
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8  Low-voltage  I-V  curves  for  Pb-Cu-7. 
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Figure  4.9  More  low-voltage  I-V  curves  for  Pb-Cu-7. 
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zero-voltage  axis  up  to  the  voltage  V-jQwer  where  the  curve  becomes  a 
straight  line.  At  lower  temperatures  there  is  a  discontinuous  jump  in  the 
voltage  at  the  critical  current,  but  the  curves  are  still  single-valued. 

At  still  lower  temperatures  the  curve  becomes  hysteretlc;  i.e..  If  the 
current  Is  Increased  from  zero  the  voltage  jumps  to  a  finite  value  at  the 
critical  current  I  ,  but  if  the  current  is  slowly  decreased  from  above  I 

C  v 

the  voltage  does  not  return  to  the  zero-voltage  axis  until  a  lower 
current  IR  is  reached.  Each  of  these  temperature  regimes  will  be  dis¬ 
cussed  in  detail. 

Single-valued,  continuous  regime 

In  this  regime  the  I-V  curve  rises  from  the  zero-voltage  axis  with 
an  initially  upward  curvature;  at  the  inflection  point  v-jnf|ection  the 
curvature  changes  sign;  when  the  voltage  Vlower  is  reached  the  I-V  curve 
becomes  a  straight  line.  The  temperature  dependence  of  these  voltages 
inflection  and  Vlower*  as  we^  as  the  temperature  dependence  of  the 
currents  Inflection  and  iower  at  which  these  volta9es  occur,  are  shown 
in  Figures  4.10  and  4.11.  It  was  found  that  within  the  accuracy  of  the 

data  the  voltages  Vinfiection  and  Vlower  are  indePendent  of  temperature. 

On  the  other  hand,  the  corresponding  currents  Inflection  and  1 lower 
seem  to  be  decreasing  with  Increasing  temperature.  (See  the  Measurement 
Techniques  chapter  for  a  discussion  on  the  measurement  of  V-|ower  and 

inflection*^  Mote  that  1f  the  currents  inflection  and  'lower  are 
multiplied  by  the  temperature-dependent  (near  Tc)  low-voltage  resistance 

(Figure  4.35)  the  resulting  voltages  are  still  decreasing  with  increasing 

temperature. 
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Figure  4.11  Temperature  dependence  of  I1ower  and  ^nf|ection  for 
Pb-Cu-13. 


Single-valued,  discontinuous  regime 

At  Intermediate  temperatures  the  I-V  curve  leaves  the  zero-voltage 
axis  dlscontinuously,  jumping  up  to  a  continuous  curve  whose  shape  Is 
very  much  like  the  shape  of  the  I-V  curve  in  the  continuous  regime.  At 
the  high-temperature  end  of  this  regime  the  I-V  curve  has  an  Inflection 
point  at  the  same  voltage  as  In  the  continuous  regime.  As  the  tempera¬ 
ture  is  lowered,  the  magnitude  of  the  voltage  jump  increases;  at  the 
low- temperature  end  of  this  regime  this  magnitude  has  Increased  beyond 

inflection  an(1  curve  no  1°n9er  has  an  inflection  point.  Just  before 
the  magnitude  cf  the  jump  reaches  Vlower  the  I-V  curve  becomes  hysteretic. 

Hysteretic  regime 

At  low  temperatures  the  I-V  curve  is  hysteretic;  this  hysteresis  can 
be  characterized  by  a  return  voltage  VR  and  a  return  current  IR.  Figures 
4.12  and  4.13  show  the  temperature  dependence  of  IR  and  VR.  It  was  found 
that  at  low  temperatures  both  IR  and  VR  become  independent  of  temperature. 
It  can  be  seen  in  Table  1  that  the  low- temperature  limit  of  VR  is  very 
nearly  equal  to  VlQwer  for  the  samples  studied.  It  was  also  found  that 
after  IR  and  VR  cease  to  depend  on  temperature  the  only  change  with 
temperature  of  the  I-V  curve  was  the  critical  current— the  finite- 
voltage  portions  of  the  I-V  curve  do  not  change  with  temperature.  Note 
that  while  the  low-voltage  resistance  Is  temperature  independent  below 
about  0.7  T/T  (Figure  4.35),  the  I-V  curve  becomes  temperature  inde- 

V 

pendent  only  below  about  0.5  T/Tc. 

The  temperature  at  which  the  return  voltage  and  the  return  current 
become  constant  is  related  to  the  sample  length;  i.e.,  longer  samples 
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Figure  4.12  Temperature  dependence  of  the  return  voltage  for  three 
Pb-Cu-Pb  samples. 
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Figure  4.13  Temperature  dependence  of  the  return  current  IR  for  the 
samples  as  in  Figure  4.12. 
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reach  this  stage  at  lower  temperatures.  In  fact,  for  some  of  the  longer 
samples  we  could  not  lower  the  temperature  low  enough  for  this  return 
voltage  to  become  constant,  although  the  trend  In  this  voltage  Indicated 
that  It  would  become  constant.  For  very  long  samples  no  hysteresis  was 
observed  in  solely  dc  measurements,  but  hysteresis  could  be  Induced  in 
these  very  long  samples  by  Irradiation  with  microwaves. 

4.2.2  Theoretical  considerations 

Hysteresis 

The  hysteresis  found  in  these  SNS  samples  is  unusual  in  that  all  of 

the  other  popular  types  of  superconducting  weak  links  (oxide  tunnel 

junctions,  point  contacts,  and  microbridges  with  superconducting  bridges) 

do  not  exhibit  temperature- independent  return  currents  at  temperatures 

where  the  critical  current  is  not  also  temperature  independent.  Thus  the 

unusual  nature  of  the  observed  hysteresis  is  due  to  the  fact  that  the 

bridge  is  made  from  a  normal  metal. 

There  have  been  three  sets  of  theories  advanced  as  explanations  for 

hysteresis  in  superconducting  weak  links.  The  first  theory  proposed  was 

that  the  capacitance  between  the  electrodes  was  the  cause  of  the 
25  26 

hysteresis.  *  This  explanation  can  be  ruled  out  for  two  reasons: 

the  planar  geometry  of  our  samples  has  far  too  little  capacitance  to  cause 

any  hysteresis,  and  these  capacitance  theories  require  a  change  in  the 

return  current  if  the  critical  current  changes. 

The  second  possible  possible  explanation  is  the  joule  heating  theory 
27 

of  Skocpol  e£  ajL  This  theory  is  also  unable  to  account  for  the 
observed  temperature  dependence  of  IR,  since  this  theory  requires  IR  to 


be  proportional  to  (1  -T/TCP;  our  samples  have  an  IR  which  Is  Independent 

of  temperature  at  low  temperatures .  In  addition.  It  will  be  shown  later 

that  the  response  of  these  microbridges  to  microwave  radiation  precludes 

the  existence  of  sufficient  joule  heating  to  cause  hysteresis. 

The  third  possibility,  and  the  one  which  seems  to  be  applicable  to 

our  samples,  is  that  the  hysteresis  is  a  result  of  an  intrinsic  relaxation 

time.  This  relaxation-time  explanation  was  first  proposed  by  Song  ,  and 

29  20 

it  has  been  put  on  a  stronger  theoretical  footing  by  other  authors.  * 
These  theories  are  set  in  a  time-dependent  Ginzburg -landau  framework,  and 
this  framework  allows  some  Qualitative  comparisons  to  be  made. 


Time-dependent  Ginzburg-Landau  relaxation  time 

The  time-dependent  Glnzburg-landau  (TDGL)  theory  contains  an  intrinsi 

relaxation  time  which  is  valid  for  superconductors  above  T  :  we  can  use 

this  time  by  considering  a  normal  metal  to  be  a  superconductor  with  a  Tc 

31 

of  0  K.  This  TDGL  relaxation  time  is  given  by 


_ nil  _  TrfS 

«■  8k(t‘  Vbrldge1  "  8kT 

where  we  have  used  s  0.  If  the  spatial -gradient  terms  of  the 

TDGL  equations  are  taken  into  account,  the  effective  relaxation  time  for 
a  microbridge  of  length  L  is  given  by 
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This  time  is  essentially  the  time  required  for  electron  pairs  to  diffuse 
through  the  normal -metal  bridge.  It  is  also  the  time  for  any  electron  to 
diffuse  through  the  normal  metal.  This  diffusion  process  is  the  reason 


why  the  effective  time  is  proportional  to  the  square  of  the  length  of  the 
bridge. 

There  are  two  important  features  of  this  relaxation  time  which  are 
easily  accessible  to  experimental  measurement.  The  first  is  that  is 
Independent  of  temperature,  and  the  second  Is  the  dependence  of  on  l. 

Baratoff  and  Kramer  (BK)  have  numerically  solved  the  TDGL  equations 
for  the  case  of  a  weak  link  whose  bridge  is  a  gapless  superconductor  with 
reduced  length  l  -  1.  They  found  that  the  I-V  characteristics  could  be 
classified  by  a  quantity  T: 

r  *  Wtj  w,,ere  tj  *  2inc  • 

C  0 

For  our  samples  t  is  a  constant,  independent  of  temperature;  thus  r 
Is  Inversely  proportional  to  I  for  our  samples. 

V* 

BK  have  shown  that  the  relaxation  time  affects  the  shape  of  the  I-V 
curve  in  two  ways.  First,  when  the  microbridge  is  in  a  finite-voltage 
state  the  magnitude  of  the  oscillating  supercurrent,  |I  |,  is  reduced 
below  the  zero-voltage  value  Ic.  For  a  constant  temperature  (correspond¬ 
ing  to  a  fixed  r  and  a  fixed  Ic)  Usl/Ic  becomes  smaller  as  Vdc  is 
Increased.  This  decrease  in  I§  as  I  is  Increased  produces  an  I-V  curve 
which  becomes  a  straight  line  at  a  smaller  value  of  I/I  than  the  value 

V# 

predicted  by  the  shunted  junction  model.  In  fact,  BK  have  shown  that 
above  the  voltage  V-jQwer  (the  dc  voltage  at  which  the  I-V  curve  becomes 
a  straight  line),  the  oscillating  supercurrent  is  almost  completely 
suppressed. 

BK  predicted  that  if  r  Increases,  V]0wer/^CR0  should  decrease.  We 
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find  this  to  be  the  case  in  our  samples,  except  perhaps  close  to  T  .  As 
the  temperature  is  lowered,  and  r  increases  since  r  is  inversely  propor¬ 
tional  to  Ic  in  our  samples,  we  find  that  VjQwer  is  constant.  (That  is, 
in  effect,  a  decrease  in  V-|0wer/IcR0  since  Ic  increases  with  decreasing 
temperature.)  This  temperature  independence  of  Vlowgr  can  be  seen  in 
Figure  4.14  and  in  Figures  4.7  through  4.9. 

The  second  relaxation  effect  found  by  BK  is  that  for  r  >  38  the  I-V 
curve  should  be  hysteretic;  this  hysteresis  is  observed  in  our  samples. 

We  have  defined  a  quantity  rR  by  the  value  of  r  at  which  our  samples  have 
just  fully  developed  hysteresis;  i.e.,  rR  «  ( 2e/fi ) I RRQ  xgff.  Values  of 
rR  for  different  samples  are  shown  in  Table  1;  although  there  is  con¬ 
siderable  variation  in  l  and  ICR0»  rR  is  found  to  be  very  weakly  depend¬ 
ent  on  length,  showing  a  slight  tendency  to  increase  as  l  decreases.  As 
was  found  for  V]ower>  it  has  been  observed  that  at  low  temperatures  the 
voltage  VR  at  which  the  I-V  returns  to  the  zero-voltage  axis  is  independent 
of  temperature.  We  have  found  that  V-|0wer  and  VR  have  nearly  the  same  low 
temperature  values. 

In  summary,  we  have  found  that  for  dc  voltages  across  the  sample 

greater  than  about  5  ft/Zer^,  the  I-V  curve  is  almost  independent  of 

temperature.  Considering  the  theoretical  work  of  BK,  we  have  concluded 

that  as  the  dc  voltage  exceeds  fi/2eT  ^  the  magnitude  of  the  oscillating 

supercurrent  begins  to  decrease;  when  Vd{.  >  5fi/2eTeff,  the  oscillating 

supercurrent  is  almost  completely  suppressed.  Further,  the  temperature 

independence  of  t  ^  is  confirmed  by  the  temperature  independence  of  VR 

2 

and  ower *  and  Teff  found  scale  nearly  as  i  . 


I  (microamps) 


Figure  4.14  I-V  curves  for  Pb-Cu-13.  In  this  figure  the  curves  for 
different  temperatures  are  not  offset  in  voltage  as  they 
are  in  previous  figures.  The  curves  were  taken  at 
temperatures  6.63,  6.34,  5.96,  5.44,  4.76,  4.34,  4.04,  and 
3.69  K. 
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Our  experimental  findings  on  hysteresis  are  not  in  complete  agreement 
with  the  theory  of  BK,  since  they  predicted  an  onset  of  hysteresis  at 
r  *  38,  while  our  result  is  that  hysteresis  beings  at  r  -  2. 

The  conclusion  that  above  Vlower  the  supercurrent  is  almost  completely 
suppressed  does  not  mean  that  for  dc  voltages  above  V-|ower  the  bridge  can 
always  be  treated  simply  as  a  normal  metal.  It  will  be  shown  in  the  next 
part  of  this  chapter  that  microwave  radiation  can  cause  a  supercurrent  to 
be  re-established. 

32 

The  theories  of  BK  and  others  also  do  not  describe  the  straight- 
line  portion  of  the  I-V  completely  satisfactorily  since  they  do  not 
predict  the  observed  behavior  of  the  "excess  current".  This  "excess 
current"  is  seen  in  Figures  4.7  through  4.9  by  the  fact  that  the  straight- 
line  portion  does  not  pass  through  the  origin;  this  portion  of  the  I-V 
is  offset  in  current  by  an  "excess  current".  These  TDGL  theories  contain 
an  "excess  current"  which  depends  on  I  ,  while  at  low  temperatures  our 

V 

samples  have  an  "excess  current"  which  is  independent  of  Ic- 

Although  the  BK  theory  does  not  agree  with  all  of  our  experimental 
findings,  the  BK  calculations  did  not  use  exactly  the  same  model  condi¬ 
tions  as  are  present  in  our  samples;  numerical  calculations  specific  to 
long  SNS  microbridges  are  needed  to  see  if  the  TDGL  equations  are  capable 
of  quantitatively  predicting  the  observed  hysteresis  and  "excess  current". 

Theories  of  the  low-voltage  curvature  of  the  I-V 

The  TDGL  theories  have  a  severe  problem  in  that  they  do  not  predict 
the  increase  in  I  with  microwave  radiation  (to  be  discussed  in  the  next 
part  of  this  chapter).  Secondly,  they  do  not  predict  the  inflection  point 
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seen  In  the  I-V  at  high  temperatures.  Theories  have  been  recently  advanced 

33  34 

to  explain  this  Initial  upward  curvature  of  the  i-v.  *  These  theories 
are  related  to  the  TDGL  theories,  but  they  Include  the  effect  of  non- 
equilibrium  electrons  in  changing  the  magnitude  of  the  gap  in  the  density 
of  states  in  the  bridge  region.  The  net  prediction  of  these  theories  is 
that  for  I  just  above  Ic  the  gap  and  hence  the  oscillating  supercurrent 
Increases  with  I,  thus  giving  rise  to  an  upward  curvature.  But  because 
our  SNS  microbridges  cannot  have  a  gap  in  the  bridge,  these  theories  are 
not  immediately  applicable  to  our  samples.  It  Is  possible  that  a 
similar  effect  could  be  found  in  the  correlation  picture  of  superconduc¬ 
tivity  which  is  applicable  to  SNS  structures.  The  relaxation  time 
used  above  is  a  time  which  governs  the  ability  of  the  correlation  function 
F  to  respond  to  temporal  changes  in  the  current,  and  the  temperature 
independence  of  this  time  may  be  the  cause  of  the  temperature  independence 
of  the  inflection  voltage  V-jnfjection'  Thls  enhanced_corre1at:'ion  explana¬ 
tion  will  be  discussed  in  more  detail  in  the  next  part  of  this  chapter, 
where  the  increase  in  Ic  by  microwave  radiation  will  be  considered. 

4.2.3  Previous  experimental  work 

There  have  been  only  three  other  measurements  of  the  low-voltage 

35 

I-V  characteristics  of  SNS  junctions.  Clarke  measured  his  junctions 
up  to  only  about  20  nV.  He  did  not  observe  any  hysteresis.  His  smaller 
junctions  had  I-V's  which  could  be  fit  by  the  shunted-junction  model, 
while  his  larger  junctions  had  features  characteristic  of  flux-flow. 

The  other  two  sets  of  authors2*’36  did  not  describe  their  I-V  measure¬ 
ments  in  detail,  and  no  hysteresis  or  unusual  curvature  was  noted. 
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4.3  Effects  of  Applied  Microwave  Radiation 

The  previous  two  parts  of  this  chapter  were  concerned  with  the  I-V 
characteristics  produced  when  the  sample  was  driven  by  only  a  dc  current. 
In  this  part  we  will  consider  the  I-V  curve  generated  when  the  sample 
was  driven  by  an  ac  current  as  well  as  a  dc  current;  ac  frequencies 
between  1  GHz  and  26  GHz  were  used.  The  most  Important  result  Is  that 
the  features  of  the  I-V  curves  can  be  classified  by  whether  or  not  the 
applied  rf  period  was  above  or  below  a  critical  period.  It  has  been 
empirically  observed  that  this  critical  period  was  determined  by  the 
Josephson  period  or  the  Ginzburg-Landau  relaxation  time,  whichever  was 
longer. 

There  were  two  basic  changes  in  the  I-V  curve  when  an  ac  current 
was  applied.  First,  the  critical  current  increased,  an  effect  which  is 
popularly  known  as  the  Wyatt-Dayem  effect.  Second,  constant-voltage 
steps  appeared  on  the  I-V  curves,  an  effect  known  as  the  ac  Josephson 
effect.  Both  of  these  effects  can  be  seen  in  Figure  4.15.  Although 
these  two  effects  are  not  entirely  unrelated,  they  will  be  presented 
separately. 

4.3.1  Previous  experimental  work  on  microwave-enhanced  critical 
currents 

The  first  observations  of  a  microwave-enhanced  critical  current 
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were  made  by  Wyatt  jit  al_.  and  by  Dayem  and  Wiegand.  They  used 

samples  of  the  uniform-thickness  constriction  type  and  observed  an 

39 

enhancement  only  very  close  to  T  .  Gregers-Hansen  et.  al_.  found  a 

V 

similar  behavior  In  scribed  tin  and  indium  bridges,  with  enhancements 

40 

only  seen  close  to  T  .  Fjordbtfge  et  al .  observed  an  enhancement  in 

V 
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Figure  4.15  I-V  curves  for  Pb-Cu-13  with  10.0  GHz  applied  microwave 
radiation.  T  =  5.44  K  and  w  -r/w,  =  1.94. 
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aluminum  point  contacts;  they  also  found  the  enhancement  to  occur  most 

strongly  very  close  to  T  . 

41-43 

Notarys  et.  al_.  conducted  a  study  of  the  dependence  of  the 
critical -current  enhancement  on  the  length  of  the  bridge.  They  used 
samples  In  which  the  weak  region  was  obtained  by  underlying  a  section  of 
a  strip  of  superconductor  with  a  lower  T  material  (gold  under  tin  or 

V 

titanium  under  tantalum)  or  by  ion-implanting  a  section  of  the  super¬ 
conducting  strip  (copper  or  iron  into  tantalum).  They  found  large  enhance¬ 
ments  in  the  tin-gold  samples,  with  the  largest  enhancements  occurring 
in  the  shortest  samples.  The  tantalum-based  samples  had  smaller  enhance¬ 
ments  than  the  tin-gold  samples,  with  the  longer  tantalum-based  samples 
having  larger  enhancements  than  similar  shorter  tantalum-based  samples. 

(The  tantalum-based  samples  had  much  longer  reduced  lengths  than  the  tin- 
gold  samples.)  For  all  of  their  samples  the  critical  current  was  enhanced 

only  at  temperatures  within  about  0.2  of  T/Tc. 

44 

Latyshev  and  Nad  conducted  a  thorough  study  of  the  frequency 
dependence  of  the  enhancement  in  variable-thickness  tin  bridges.  They 
found  a  lower  frequency  limit  which  was  essentially  independent  of 
temperature.  They  also  found  an  upper  limit  for  the  frequency;  if  the 
frequency  of  the  applied  microwaves  was  higher  than  the  gap  frequency, 
they  found  no  enhancement. 

It  will  be  shown  below  that  the  results  of  all  of  these  previous 
experiments  differ  significantly  from  our  results. 

4.3.2  Microwave-enhanced  critical  current,  experimental  results 

Figures  4.16  through  4.19  show  the  temperature  dependence  of  the 
critical  current  Ic  when  a  microwave  field  was  applied  for  several  samples. 
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Figure  4.16  ICR0  versus  T/Tc  for  Pb-Cu-13  with  three  different 

frequencies  of  applied  microwave  radiation.  In  each  case 
the  data  were  taken  at  the  constant  rf  power  level  which 
maximized  IrRg  at  low  T.  The  solid  line  Is  the  theoreti¬ 
cally  calculated  value  for  i  »  8.9,  but  with  magnitude 
reduced  by  a  factor  of  0.7. 
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Figure  4.17  IcRg  versus  T/Tp  for  Pb-Cu-B2  with  five  different  fre¬ 
quencies  of  applied  microwave  radiation.  For  each 
frequency  the  rf  power  level  was  held  constant.  The 
solid  line  Is  the  theoretical  prediction,  reduced  by  a 
factor  of  0.7,  for  l  -  16.5. 


-  i 
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Figure  4.18  IcRo  versus  T/Tc  for  Pb-Cu-B3  with  10  GHz  radiation;  the 
data  was  taken  at  a  constant  rf  power  level.  The  solid 
line  is  the  calculation  for  %  *  30  which  is  the  estimated 
value  for  this  junction.  No  critical  current  could  be 
observed  In  the  absence  of  microwaves. 


versus  T/Tc  for  Pb-Cu-B2  with  10  GHz  radiation  of 
erent  power  levels.  18  dB  is  the  power  level  that 
produced  the  first  maximum  in  ICR0;  7.5  dB  produced  the 
second  maximum. 


Figure  4.19  ICR0 
dlff 


It  can  be  seen  that  the  relative  enhancement  I  /I  (I  is  the  critical 

V  w  U 

current  in  the  absence  of  microwave  radiation)  depended  on  the  temperature, 
the  frequency,  and  the  length  of  the  bridge.  We  have  found  that  there 
was  a  minimum  frequency  necessary  to  observed  an  enhanced  I  .  Two 
empirical  criteria  have  been  established  for  this  minimum  frequency;  it 
was  necessary  to  satisfy  only  one  of  the  criteria  to  observe  an  enhance¬ 
ment. 

The  first  criterion  was  that  had  to  be  greater  than  the  Oosephson 
frequency  u>j  3  2eIcR0/fi.  For  a  given  microwave  frequency  this  determined 
a  temperature  above  which  an  enhancement  could  be  seen,  provided  that  the 
enhanced  Ic  was  large  enough  to  be  observed  above  the  noise. 

The  second  criterion  was  that  >  A/t^  (A  %  2),  where  is  the 
effective  time-dependent  Ginzburg-Landau  relaxation  time  discussed  in 
detail  in  the  second  part  of  this  chapter.  It  is  important  to  note  that 
Teff  indePendent  temperature;  thus  when  u>rf  was  significantly  greater 
than  A/Teff,  Ic  enhancement  would  occur  at  all  temperatures  regardless 
of  whether  or  not  u>rf  >  u>j.  Since  Te^  is  proportional  to  l  U  is  the 
reduced  length  of  the  bridge),  this  second  criterion  was  satisfied  for 
longer  bridges  at  lower  frequencies.  It  can  also  be  observed  that  as 
became  much  greater  than  u>j,  or  as  became  much  greater  than 
A/te^»  the  relative  enhancement  became  much  greater,  becoming  as  much 
as  two  orders  of  magnitude  for  the  longest  bridges  studied. 

Figure  4.19  Illustrates  the  effect  that  different  power  levels  of 
a  constant  microwave  frequency  had  on  the  critical  current.  It  can  be 
seen  that  the  effect  of  the  microwave  field  can  be  interpreted  as  having 
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changed  the  temperature  dependence  of  I  R.  As  the  power  was  increased,  the 

v 

ICR  versus  T  curves  took  on  shapes  appropriate  for  smaller  values  of  i. 

This  slope-changing  effect  can  also  be  seen  in  Figures  4.16  and  4.17.  The 
longer  the  l,  or  the  higher  the  frequency,  the  more  pronounced  the  change. 
However,  while  the  microwave  radiation  radically  changed  the  slope  of  the 
ICR  curves,  the  low-temperature  values  of  ICR  remained  at  low  levels.  At 
no  time  did  we  observe  an  enhanced  I  R  larger  than  the  T=0  value  predicted 
for  a  given  junction  by  the  Usadel  equations. 

Our  experimental  results  are  quite  different  from  those  of  previous 
workers.  We  find  that  the  critical  current  can  be  enhanced  at  all 
temperatures  below  the  transition  temperature  of  the  banks,  provided 
that  the  applied  microwaves  have  a  high  enough  frequency.  This  disagree¬ 
ment  Is  probably  due  to  the  fact  that  our  samples  had  bridges  composed 
of  a  material  which  Is  not  superconducting  at  any  temperature,  whereas 
the  previous  workers  had  bridges  which  were  superconducting  at  relatively 
high  temperatures. 

4.3.3  Theoretical  considerations  on  microwave-enhanced  critical 
currents 

There  have  been  two  different  theories  advanced  to  explain  the 

microwave-enhancement  of  the  supercurrent  in  superconducting  microbridges. 

Both  of  these  theories  are  based  on  the  effect  that  non-paired  electrons 

(quasi particles)  have  on  the  density  of  paired  electrons.  The  first 

45 

theory  was  originally  advanced  by  Eliashberg  and  was  further  developed 

46 

by  Chang  and  Scalaplno.  This  theory  suggests  that  the  microwaves 
change  the  quasiparticle  distribution  by  Increasing  the  average  energy 
of  the  quasi particles;  this  change  results  in  an  Increase  in  the 
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superconducting  energy  gap.  Since  the  number  of  quasi particles  is  highest 

near  T  ,  the  gap  is  increased  most  near  T  .  It  seems  that  this  mechanism 
»»  c 

should  not  apply  to  our  samples,  since  there  is  no  gap  in  the  normal  metal 
bridge.  An  increase  in  the  gap  of  the  banks  could  be  induced,  but  since 
Ic  depends  at  most  on  (^ap^)  »  the  gap  in  the  banks  would  have  to  be 
increasing  by  orders  of  magnitude  at  temperatures  well  below  T  .  This 

V 

theory  predicts  only  a  very  small,  if  any,  change  in  the  gap  at  such 
temperatures.  Thus  it  appears  that  this  gap-enhancement  theory  does  not 
apply  to  our  samples. 

The  second  theory  was  first  proposed  by  Lindelof^  and  has  been 

40  . 

developed  on  a  more  theoretical  basis  by  Aslamazov  and  Larkin.  This 
theory  supposes  that  since  quasi  particles  of  energy  less  than  the  gap 
are  trapped  in  the  bridge,  the  effect  of  the  microwaves  is  to  broaden 
the  spatial  distribution  of  the  quas i parti  cl es .  Since  charge-neutrality 
must  be  preserved,  the  spatial  distribution  of  paired  electrons  is  also 
broadened.  This  idea  is  depicted  schematically  in  Figure  4.20.  This 
spatial  broadening  results  in  an  increased  pair  density  in  the  center  of 
the  bridge,  thus  Increasing  the  critical  current. 


The  idea  that  it  is  the  change  in  the  quasi  particles  which  is  causing 
the  increased  critical  current  is  supported  by  our  observation  of  two 
criteria  for  enhanced  superconductivity.  If  the  microwave  period  is  longer 
than  both  the  Josephson  period  and  the  relaxation  time  r^.,  the  ac 
current  Is  carried  almost  completely  by  the  paired  electrons;  the  quasi¬ 
particles  are  essentially  unaffected.  But  if  the  ac  current  oscillates 
faster  than  either  of  these  times,  the  quasi particles  are  forced  to  carry 
some  of  this  ac  current  and  thus  their  spatial  distribution  is  altered. 
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This  second  theory  does  not  predict  an  enhanced  I  in  all -super- 
conducting  mi crobri dges  far  from  T  .  In  this  case,  the  pair  density 
is  flat  in  the  center  of  the  bridge,  and  spatially  broadening  the  pair 
density  would  not  result  in  an  increased  pair  density  in  the  center  of 
the  bridge.  Thus,  this  theory  can  explain  why  our  samples  show  enhanced 
critical  currents  at  all  temperatures,  while  all -superconducting  micro- 
bridges  only  show  an  enhanced  I  near  T  . 

The  gap-enhancement  theories  also  predict  an  enhancement  only  near 
Tc  for  all -superconducting  microbridges.  This  is  due  to  the  fact  that, 
in  contrast  to  our  SNS  samples,  the  relaxation  time  for  all -superconduct¬ 
ing  microbridges  gets  shorter  as  the  temperature  is  lowered.  The  gap- 
enhancement  theories,  however,  have  not  been  very  successful  in  explain¬ 
ing  the  fact  that,  at  least  in  some  of  the  all -superconductor  systems 
studied,  longer  bridges  show  a  smaller  enhancement  than  shorter  bridges. 

This  spatial -averaging  theory  might  also  account  for  the  low-voltage 
curvature  in  the  dc  I-V  characteristics.  In  this  region  of  the  I-V, 
there  is  an  oscillating  voltage  across  the  bridge,  and  thus  there  is  an 
oscillating  quasiparticle  current.  This  oscillating  quasiparticle 
current  might  broaden  the  spatial  distribution  of  the  paired  electrons 
In  the  bridge,  thus  increasing  the  effective  critical  current.  As  the 
dc  voltage  across  the  bridge  is  Increased,  the  quasiparticles  oscillate 
faster  and  the  critical  current  is  increased.  This  increase  in  Ic  as 
the  voltage  is  raised  causes  an  upward  curvature  in  the  I-V.  This 
effect  only  applies  for  voltages  less  than  since  at  voltages 

above  this  value  the  oscillating  supercurrent  decreases  due  to  the 
Inability  of  the  supercurrent  to  oscillate  with  a  period  much  shorter 
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than  Tgff. 

The  major  difficulty  in  applying  these  theories  to  our  samples  is 
that  the  theoretical  treatments  have  been  done  only  for  all -superconducting 
systems.  It  is  certainly  possible  that  a  rigorous  calculation  for  SNS 
systems  might  provide  yet  another  explanation. 

4.3.4  AC  Josephs on  effect 

It  was  mentioned  in  the  introduction  that  when  there  is  a  voltage 
across  a  Josephson  junction,  the  supercurrent  oscillates  in  time.  If 
both  a  dc  and  an  ac  current  are  used  to  bias  the  junction,  the  interaction 
between  the  oscillating  supercurrent  and  the  ac  drive  current  produces 
constant-voltage  steps  in  the  dc  I-V  characteristics.  These  steps  can  be 
seen  in  Figure  4.15.  If  the  junction  is  ideal,  i.e.,  the  current  is 
sinusoidally  dependent  on  the  phase  difference  and  there  are  no  relaxation¬ 
time  effects,  the  steps  should  occur  only  at  dc  voltages  which  are  multiples 
of  the  voltage  V^  =  (fi/2e)wrf,  where  is  the  frequency  of  the  applied  ac 
current.  Furthermore,  the  step-height  (the  dc-current  length  of  the  step) 
oscillates  as  the  magnitude  of  the  ac  current  increases.  Russer  has 
calculated  on  an  analog  computer  the  step-height  versus  the  ac  current 
level  for  an  ideal  junction  for  several  values  of  wrf/wj  (“j  =  (2e/fi)IcR0) . 
Figure  4.21  shows  his  results. 

But  the  junctions  produced  in  this  work  were  not  ideal.  Even  though 
the  current-phase  relation  was  sinusoidal,  it  has  been  shown  in  the 
previous  sections  of  this  chapter  that  there  are  strong  relaxation-time 
effects  present.  The  influence  of  this  relaxation  time  on  the  I-V  curves 
when  the  sample  is  driven  with  an  ac  and  a  dc  current  will  be  examined 
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Figure  4.21  Theoretical  calculation  of  the  dependence  of  the  critical 
current  (n=0)  and  the  first  step  (n=l)  on  the  rf  current 
level  for  an  ideal  resistively-shunted  junction;  y  = 

(from  Ref.  49).  rT  J 
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for  Pb-Cu-13;  the  rest  of  the  samples  studied  showed  similar  effects. 

Sample  Pb-Cu-13  had  an  effective  relaxation  time  t  ^  of  2.8  x  10*^ 
seconds.  This  corresponds  to  a  frequency  of  5.7  GHz.  We  have  found  that 
if  the  applied  ac  current  was  much  less  than  5.7  GHz,  the  data  was  some¬ 
what  in  agreement  with  the  ideal -junction  calculations.  If  the  applied 
ac  current  was  much  higher  than  5.7  GHz,  the  data  was  quite  different 
than  the  ideal -junction  calculations. 

Figure  4.22  shows  the  I-V  curves  for  various  ac-current  levels  for 
an  applied  ac  frequency  of  1.94  GHz  (wrf/wj  =1.41).  It  can  be  seen  that 
the  steps  occur  only  at  multiples  of  (fi/2e)u)  ^  =  4.0  microvolts.  Figure 
4.23  shows  a  plot  of  the  step-height  versus  ac  current  level  for  this 
I-V  curve.  Since  Russer  had  not  calculated  the  step-heights  for 
Urf/wj  =  1.4,  we  did  a  numerical  calculation  in  the  ideal -junction  limit. 
(See  the  Appendix  for  the  details  of  this  calculation.)  Note  that  the 
theory  and  data  have  similar  curves;  in  particular,  the  magnitudes  of  the 
step  heights  agree.  But  the  period  of  oscillation  of  the  step  heights 
with  ac  current  level  do  not  agree.  This  disagreement  could  be  caused  by 
the  fact  that  the  ac  current  had  a  magnitude  such  that  regions  of  the  I-V 
curve  at  voltages  greater  than  ii/2e  t  ^  were  being  sampled.  Figure  4.24 
shows  1.94  GHz  step-height  data  for  the  same  sample,  but  at  a  higher 
temperature.  Again,  the  qualitative  behavior  is  similar  to  that  pre¬ 
dicted  by  the  ideal -junction  calculations.  In  addition,  the  I-V  curve 
at  this  temperature  only  showed  steps  at  voltages  which  were  multiples 
of  Cfi/2e)u 

Figure  4.15  shows  I-V  curves  for  10.0  GHz  applied  ac  currents.  Note 
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I-V  curves  for  Pb-Cu-13  with  1.94  GHz  applied  microwave 
radiation;  T  *  6.34  K  and  uj^/uj  *  1.41. 
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Figure  4.23  Step-height  versus  rf  current  level  for  the  critical 
current  (O'th  step)  and  the  first  step  of  the  data  in 
Figure  4.22.  The  open  circles  are  the  results  of  a  numerical 
calculation  for  an  ideal  junction  with  u>rf/u)j  *  1.41  (see 
the  Appendix  for  details).  The  data  has  been  normalized  to 
the  calculation  by  requiring  the  critical  currents  at 
I  f  =  0  to  agree  and  by  requiring  the  first  zeros  of  the  O'th 
step  to  occur  at  the  same  value  of  Irf. 


that  in  addition  to  steps  appearing  at  multiples  of  V  *  (fi/2e)wrf  (20 
microvolts),  there  are  constant-voltage  steps  at  submultiples  of  this 
voltage.  These  are  often  called  subharmonic  steps.  Figure  4.25  shows 
the  step  height  versus  ac  current  level  for  these  I-V‘s.  There  are  two 
features  which  are  not  in  agreement  with  the  predictions  of  the  analog 
and  numerical  calculations:  First,  the  variation  of  step  height  with  ac 
current  level  is  somewhat  erratic.  Second,  there  is  an  initial  increase 
in  the  critical  current,  and  the  maximum  height  of  the  first  step  exceeds 
the  zero-power  critical  current.  This  is  the  enhanced  critical  current 
discussed  previously  in  this  chapter. 

Figure  4.26  shows  I-V  curves  for  ac  currents  at  24  GHz  ( (fi/2e)a)rf  = 
50  microvolts),  and  Figure  4.27  shows  the  corresponding  step  height 
versus  ac  current  level  data.  At  this  frequency,  there  are  strong  sub- 
harmonic  steps,  and  the  step  height  versus  ac  current  level  data  is 
quite  different  from  that  expected  from  the  calculations.  In  addition, 
the  critical -current  has  been  enhanced  so  much  that  the  I-V  curve  becomes 
hysteretic. 

It  has  been  found  that  this  behavior  occurs  for  all  the  samples 
measured.  If  «  l/te^,  the  resulting  I-V  curves  with  an  ac  current 
drive  are  quite  similar  to  those  predicted  by  both  analog  and  numerical 
calculations.  If  u)r^  »  1/Tfi^,  there  are  constant-voltage  steps  at 
submultiplies  of  (-fi/2e)w  and  the  step  height  versus  ac  current  level 
data  Is  quite  different  from  that  predicted  by  the  calculations. 

It  has  also  been  observed  that,  for  Pb-Cu-13  at  least,  the  constant- 
voltage  steps  can  be  seen  at  dc  voltages  as  high  as  two  millivolts.  In 
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Irf  (arb.  units) 


Figure  4.25  Step-height  versus  rf  current  level  for  the  O'th  and 
1  * st  steps  of  Figure  4.15. 
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Figure  4.27  Step-height  versus  rf  current  level  for  the  O'th  and  1'st 
steps  of  Figure  4.26. 
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previous  discussions  in  this  chapter  it  was  speculated  that  above  V-jower 
(about  fifty  microvolts)  the  bridge  acted  as  if  it  were  completely  normal. 
It  would  thus  seem  that  there  is  a  contradiction.  In  fact,  however,  there 
is  no  contradiction,  as  the  following  argument  will  show.  Figure  4.28 
shows  the  results  of  an  ideal -junction  numerical  calculation  for 
°Vf  *  WJ*  an<*  ^or  a  vo^9e  ^  Ic*0  88  TO  (fi/2e)wj .  It  can  be  seen 
that  the  supercurrent  has  two  frequencies  of  oscillation,  one  much  faster 
than  ojj  and  one  much  slower  than  uij.  (This  is  just  a  result  of  the 
oscillating  supercurrent  beating  against  the  applied  ac  current.)  Now 
if  a  relaxation  time  equal  to  1/wj  is  introduced,  as  seen  in  Figure 
4.29,  there  are  again  two  frequencies  of  oscillation;  the  high-frequency 
oscillation  of  the  supercurrent  is  significantly  suppressed,  but  at 
some  times  there  is  still  a  significant  supercurrent.  It  is  this 
remaining  supercurrent  which  is  responsible  for  the  observation  of  the 
constant  voltage  step  at  voltages  much  higher  than 'fi/2exe^f.  (A  relaxa¬ 
tion  time  was  introduced  into  the  calculation  by  suppressing  the  oscilla¬ 
ting  supercurrent  by  a  factor  dependent  on  the  instantaneous  voltage 
across  the  bridge,  i.e., 

‘s  *  +  v2/vreux> 

where  V  is  the  instantaneous  voltage  across  the  bridge  and 

Vrelax  *  ^^eTrelax‘  See  the  APPenc^x  f°r  details  of  this  calculation.) 

It  is  important  to  realize  that  It  is  the  presence  of  the  applied 
ac  current  which  gives  rise  to  an  ac  supercurrent.  If  the  applied  ac 
current  were  absent,  there  would  be  no  supercurrent.  Thus,  the  observa¬ 
tion  of  a  constant-voltage  step  at  a  high  dc  voltage  V  does  not  imply 
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time  (units  of  Tj) 


Figure  4.28  Numerical  simulation  of  the  supercurrent  (Is)  versus  time 
for  an  ideal  resistively-shunted  junction  with  applied  rf 
radiation.  Iqc  *  10  Ic,  Irf  =  10  Ic,  ujrf  =  «j  =  ( 2e/+i )  ICR0 . 
The  time  is  measured  in'  units  of  tj  =  1/Wj. 
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Figure  4.29  Numerical  simulation  of  the  supercurrent  (I$)  versus  time 
for  a  resistively-shunted  junction  with  an  intrinsic  relaxa 
tion  time  Teff  and  an  applied  rf  current.  IDc  =  10  Ic, 

Irf  ■  10  Ic»  «rf  *  «J»  and  *  1/wj.  (See  the  Appendix 
for  the  details  of  this  simulation.) 
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that  there  Is  a  supercurrent  oscillating  with  frequency  2eV/f» . 
to  establish  that  a  supercurrent  is  oscillating  at  2eV/fi ,  it  is  necessary 
to  see  a  constant- voltage  step  with  an  applied  ac  current  of  frequency 
u>rf  *  2eV/fi. 

The  observation  of  a  high-voltage  step  does  imply  that  joule 
heating  has  not  heated  the  bridge  to  the  point  where  there  can  no  longer 
be  a  supercurrent,  since  no  amount  of  mixing  with  an  applied  ac  current 
could  then  result  in  a  reappearance  of  a  supercurrent.  Thus  the  observa¬ 
tion  of  a  step  at  2  mV  implies  that  the  hysteresis  in  the  I-V  seen  at 
50  microvolts  is  not  caused  by  heating. 

4.3.5  Previous  experimental  observations  of  the  ac  Josephson 
effect  i n  SNS  s  tructures 

35 

Clarke  applied  very  low  frequency  (<  2  MHz)  radiation  to  his 
junctions  and  observed  both  harmonic  and  subharmonic  steps.  He  attributed 
the  subharmonic  steps  to  the  large  capacitance  of  his  junctions.  This  is 
probably  correct  since  his  layered  structures  had  significant  capacitance. 
He  also  found  that  the  step-height  versus  ac  current  level  was  similar 

to  that  predicted  by  calculations. 

36 

Komarovskikh  et  al_.  observed  steps  at  39.8  GHz  in  their  planar 
junctions,  but  very  little  additional  information  was  reported. 

4.4  Dynamic  Resistances  and  the  High-Voltage  I-V 

The  three  previous  sections  of  this  chapter  dealt  with  phenomena 
associated  with  the  existence  of  supercurrents  across  the  SNS  sample. 

The  experimental  I-V  curves  also  showed  features  which  were  not  due  to 
supercurrents ,  but  which  Instead  were  due  to  the  influence  the  SNS 
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Interface  had  on  the  density  of  states  (and  the  occupation  of  these  states) 
in  the  adjoining  superconductor  and  normal  metal.  These  non-supercurrent- 
related  effects  are  the  subject  of  this  section. 

4.4.1  Experimental  results 

Figure  4.30  shows  a  dc  I-V  curve  for  Pb-Cu-13  with  the  voltage 
ranging  from  zero  to  six  millivolts.  The  supercurrent-related  features 
are  only  present  below  one  millivolt;  the  remainder  of  the  I-V  shows  only 
a  very  weak  strucutre.  Figures  4.31  and  4.32  show  dV/dl  data  for  this 
sample;  it  can  be  seen  that  the  weak  features  of  the  I-V  curve  are  very 
pronounced  on  the  dV/dl  curves.  Figure  4.33  shows  a  similar  dV/dl  curve 
for  Pb-Cu-7. 

Figure  4.34  shows  a  series  of  dV/dl  curves  for  Pb-Cu-13,  thus 
illustrating  the  temperature  dependence  of  the  various  features.  The 
7.326  K  curve  was  taken  above  the  Tc  of  the  banks;  note  that  in  all  of 
the  other  curves  the  magnitude  of  dV/dl  never  exceeds  the  normal  state 
resistance  obtained  at  7.326  K. 

There  are  two  features  of  the  dV/dl  curves  which  will  be  examined 

in  detail.  The  first  feature  is  the  flat  region  occurring  between  about 

0.1  mV  and  1  mV.  This  is  the  constant-resistance  region  between  V, _ 

i  ower 

and  Vupper  which  has  been  discussed  previously;  the  magnitude  of  this 
voltage-independent  resistance  was  previously  designated  by  RQ.  Figure 
4.35  shows  the  temperature  dependence  of  Rg.  At  low  temperatures  Rg  is 
nearly  constant,  while  close  to  Tc  the  resistance  Rg  increases  rapidly 
with  temperature. 

The  second  Interesting  feature  is  the  large  peak  in  the  dV/dl  curve 


V  (millivolts) 
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I  (milliamps) 


Figure  4.30  High-voltage  I-V  curve  for  Pb-Cu-13  at  T  *  1.996  K. 
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Figure  4.33  dV/dl  versus  I  for  Pb-Cu-7,  T  *  1.996  K. 


dV/dl  (arb.  units) 


Figure  4.34  dV/dl  versus  I  for  Pb-Cu-13  for  several  temperatures. 

The  tic  marks  on  the  right-hand  sides  of  the  figure  indicate 
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which  occurs  near  2.7  mV  at  low  temperatures.  Figure  4.36  shows  the 
temperature  dependence  of  the  voltage  at  which  this  peak  occurs.  At 
temperatures  very  close  to  T  there  is  a  second,  smaller  peak  occurring 

V 

at  about  twice  the  voltage  of  the  large  peak. 

In  addition  to  these  two  features,  there  are  smaller  bumps  occurring 
between  the  flat  region  of  RQ  and  the  large  peak.  It  will  be  argued 
below  that  these  bumps  are  probably  related  to  the  large  peak. 

4.4.2  Theoretical  considerations  and  previous  experimental  work 
The  temperature  dependence  of  Ro 

The  resistance  of  a  superconducting-normal  interface  was  first 
measured  by  Pippard  and  co-workers.50,51  They  used  SNS  structures  with 
copper  as  the  normal  metal  and  lead  or  tin  as  the  superconductor.  Their 
experimental  results  were  very  similar  to  our  measurements  of  the  tempera¬ 
ture  dependence  of  RQ. 

Pippard  et  aK  were  also  the  first  to  provide  an  explanation  for 

the  observed  data.  A  brief  summary  of  their  theory  will  be  given  here. 

52 

Andreev  first  developed  the  nature  of  the  transport  of  electrons 
across  an  SN  interface.  He  found  that  if  an  electron  in  the  normal  metal 
encounters  the  interface,  what  happens  to  it  depends  on  its  energy.  If 
the  electron  has  an  energy  much  greater  than  the  gap  of  the  super¬ 
conductor,  it  simply  travels  across  the  interface.  But  If  the  electron 
has  an  energy  less  than  the  gap  of  the  superconductor,  the  electron  does 
not  simply  cross  the  interface;  Instead,  the  electron  from  the  normal 
metal  combines  with  a  quasiparticle  in  the  superconductor  to  form  a 
condensed  pair,  and  a  hole  is  ejected  back  into  the  normal  metal.  Thus, 
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an  electron  with  energy  less  than  the  gap  is  reflected  as  a  hole  going  in 
the  opposite  direction. 

Pippard  et  al_.  pointed  out  that  If  an  electron  in  the  normal  metal 
penetrates  into  the  superconductor  it  will  be  out  of  equilibrium;  it  will 
scatter  inelastically  for  a  while,  and  eventually  combine  with  another 
quasiparticle  to  form  a  pair.  This  scattering  and  subsequent  recombina¬ 
tion  causes  energy  dissipation  in  the  superconductor,  thus  giving  rise  to 
an  extra  resistance.  In  addition,  there  will  now  be  an  electric  field  in 
the  superconductor  due  to  the  non-equilibrium  quasiparticle  distribution. 
(Artemenko  et  a]_.  have  calculated  that  this  electric  field  consists  of 
an  initial  jump  in  the  electric  field  at  the  Interface,  followed  by  an 
exponentially  decaying  electric  field  with  a  decay  length  Ab=*  (4DreT/irAP 

where  D  Is  the  diffusion  constant  and  rQ  is  the  inelastic  scattering 
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time.  Dolan  and  Jackel  directly  measured  this  exponentially  decaying 
electric  field  and  found,  for  tin,  that  A^CT)  *  0-5)0  -T/Tc)®*^ 
microns). 

The  magnitude  of  this  extra  resistance  depends  on  the  number  of 
electrons  penetrating  into  the  superconductor.  At  low  temperatures, 
nearly  all  of  the  electrons  in  the  normal  metal  will  have  energy  less  than 
the  gap,  and  there  will  be  very  little  extra  dissipation.  As  the  tempera¬ 
ture  Is  raised,  the  gap  decreases  and  the  energy  spread  of  the  electrons 
in  the  normal  metal  increases,  thus  the  manber  of  electrons  with  energy 
greater  than  the  gap  Increases.  The  number  of  electrons  penetrating  the 
superconductor  Increases,  and  the  extra  dissipation  Increases. 

Pippard  et  aj.. ,  using  a  phenomenological  model,  calculated  that 
the  resistance  of  an  SNS  junction  should  be  given  by 
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R0/R0K  3  (1+C)/£1+C  tanh(pAao/2kBT)] 

where  p  Is  a  constant,  RQK  is  the  resistance  at  T=0,  is  the  gap  in 
the  superconductor  far  from  the  interface,  and  C  *  rq(Tc)/RqK  - 1 . 

Pippard  et  al_.  found  from  their  experimental  data  that  p  was  between  1.4 
and  1.7.  Figure  4.35  shows  a  fit  of  our  data  to  this  equation,  and  the 
fit  is  good  using  p  *  1.6. 

Since  the  work  of  Pippard  et  al..,  many  authors  have  attempted  to 
improve  on  the  somewhat  phenomenological  theory  of  Pippard  et  al_.  by 
performing  rigorous,  microscopic  calculations.  All  of  these  authors 
have  computed  the  resistance  to  be  Rg  =  RQK  +  R',  where  R'  is  the  extra 

resistance.  Schmid  and  Schon5^  found  that  R'  'v  (1  -T/T  )“®*^  (for 

V 

lead  electrodes),  Artememko  et  al..55  found  R'  'v  (1  -T/Tc)"0'25,  and 
KrahenbUhl  and  Watts-Tobin56  found  R'  *  (1  -T/Tc)~0,63.  All  of  these 
microscopic  theories  should  be  applicable  for  0.85  <  T/Tc  <  .99. 

We  have  fit  our  data  to  the  equation  r0~rok+r-|0  -T/Tc)~p.  We 
have  found  that  the  best  fit  is  for  p  =0.11,  as  seen  in  Figure  4.36. 

The  difficulty  with  this  fit  is  that  this  gives  a  value  of  RgK  which  is 

negative,  a  result  which  is  clearly  not  physically  possible.  We  thus 
conclude  that  our  data  cannot  be  fit  by  this  equation.  The  reason  for 
the  inability  of  the  microscopic  theories  to  describe  our  data  Is  not 
known.  It  is  possible  that  the  microscopic  theories  do  not  apply  because 
they  were  worked  out  for  a  single  SN  interface;  they  do  not  include  the 
effect  that  one  of  the  superconducting  banks  has  on  the  other,  as  a 
result  of  the  change  in  the  occupied  density  of  states  of  the  normal  metal. 
(The  phenomenological  model  of  Pippard  et  al_.  did  contain  a  two-interface 
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SNS  junction,  however.) 

High-voltage  structure 

The  peak  in  the  dV/dl  curve  which  occurs  at  high  voltage  (V  %  2  mV) 
was  not  observed  in  any  of  the  previous  studies  of  SNS  structures,  since 
these  studies  did  not  measure  the  I-V  at  such  high  voltages.  However,  a 
similar  peak  in  the  dV/dl  curve  of  SN  point  contacts  has  been  observed 
by  Kaiser-Dieckhoff.^ 

CO 

Wang  et  al_.  observed  peaks  in  the  dV/dl  curves  of  niobium  micro¬ 
bridges  at  a  voltage  corresponding  to  the  gap  in  the  bridge.  They 
attributed  this  peak  to  self-coupling  to  the  Josephson  oscillation  at 

high  temperatures  and  tothe  presence  of  a  hot-spot  in  the  bridge  at  lower 
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temperatures.  Gregers-Hansen  et  al_.  observed  peaks  in  the  dV/dl  curves 
of  planar  indium  microbridges  at  voltages  corresponding  to  subharmonics 
of  the  gap  voltage.  They  attributed  this  structure  to  pair-breaking  by 
subharmonics  of  the  Josephson  oscillation.  The  peak  in  our  dV/dl  curves 
cannot  be  explained  by  coupling  to  Josephson  radiation  since  at  these 
high  voltages  there  was  no  oscillating  supercurrent  observed  in  our 
microwave  experiments.  The  possibility  that  heating  is  the  cause  of  our 
peak  will  be  discussed  below. 

Since  none  of  the  theories  of  SNS  junctions  or  SN  interfaces  have 
considered  such  high  voltages,  we  cannot  look  to  those  theories  for  an 
explanation.  Instead  a  somewhat  speculative  explanation  will  be  given 
here. 

It  has  been  mentioned  above  that  electrons  from  the  normal  metal 
which  penetrate  Into  the  superconductor  cause  an  extra  resistance,  and 
the  magnitude  of  this  extra  resistance  depends  on  the  number  of  electrons 
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In  the  normal  metal  which  have  an  energy  greater  than  the  energy  gap  in 
the  superconductor.  Now  if  a  voltage  V  is  established  across  the  bridge, 
the  mean  electron  energy  (with  respect  to  the  fermi  energy)  at  the  low- 
energy  end  of  the  bridge  will  be  higher  than  it  would  be  if  V  *  0.  (This 
is  due  to  the  fact  that  for  our  copper  bridges  the  mean  free  path  was 
about  15  nm;  since  the  inelastic  collision  time  for  copper  at  low  tempera¬ 
tures  is  about  IQ'10  seconds®®,  this  gives  a  mean  diffusion  distance  between 
inelastic  colli  sons  of 


This  distance  is  longer  than  the  length  of  our  normal  metal  bridges.) 

Thus,  as  the  voltage  across  the  bridge  increases,  the  number  of  electrons 
penetrating  into  the  superconductor  increases,  as  illustrated  in  Figure 
4.38.  In  Figure  4.38a,  most  of  the  electrons  in  the  normal  metal  bridge 
are  below  the  gap,  and  the  extra  resistance  is  low.  In  Figure  4.38b, 
the  peak  in  the  electron  distribution  is  nearly  at  the  gap,  and  the  extra 
resistance  is  increased.  In  Figure  4.38c,  the  peak  in  the  energy  distri¬ 
bution  is  far  above  the  gap.  (In  fact,  the  number  of  electrons  with 
energy  just  above  the  gap  is  less  in  (c)  than  in  (b).)  But  now,  since 
most  of  the  electrons  are  entering  the  superconductor  far  above  the  gap, 

the  extra  resistance  is  decreased  from  (b)  because  high-energy  quasi- 

64 

particles  combine  faster  than  low-energy  quasi  particles.  (This  increase 

in  the  recombination  rate  is  due  to  the  fact  that  the  phase  space 

2 

available  to  the  emitted  phonon  varies  as  fl  ,  where  SI  is  the  phonon 
energy.  If  a  quasi  particle  of  energy  A  combines  with  another  quasi - 
particle  with  energy  A,  the  available  phase  space  is  proportional  to 


Voltage  of  Peak  in  dV/dl  (mV) 
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Figure 


T/Tq 


.37  Temperature  dependence  of  the  voltage  at  which  the  peak 
In  dV/dl  occurs  for  Pb-Cu-13.  The  solid  line  is  the  gap 
voltage  for  Pb  divided  by  2.73.  The  data  has  been  divided 
by  2.825. 
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(a) 

V=0 


(b) 
V  =  A 


(c) 

Vs  2  A 


Ef 


Figure  4.38  Schematic  representation  of  the  spatial  variation  of  the 
energy  of  the  peak  in  the  electron  distribution  across  a 
SNS  microbridge.  A  is  the  gap  in  the  superconductor,  Ef 
is  the  fermi  energy,  and  Vis  the  voltage  across  the  micro' 
bridge. 
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2  2 

(2 A)  a  4A  .  If  one  of  the  quasi  particles  has  energy  2A  and  the  other  A, 

2  2 

the  available  phase  space  is  proportional  to  (3A)  *  9A  .  Thus  a  quasi  - 
particle  with  energy  2A  combines  more  than  twice  as  fast  as  one  with 
energy  A.) 

Thus,  as  the  voltage  increases,  the  extra  resistance  increases  and 
then  decreases.  Figure ‘4.37  shows  the  temperature  dependence  of  the 
voltage  of  the  dV/dl  peak  along  with  the  temperature  dependence  of  the 
energy  gap.  The  shape  is  in  fair  agreement,  but  the  voltage  of  the  peak 
is  larger  than  the  gap  voltage.  The  reason  for  this  discrepancy  is 
possibly  the  following:  the  peak  In  dV/dl  should  occur,  according  to 
the  above  model,  when  the  energy  of  the  peak  in  the  electron  distribution 
at  the  SN  interface  Is  equal  to  the  gap  in  the  superconducting  banks. 

For  this  to  happen,  the  voltage  drop  across  the  normal  metal  bridge  must 
be  equal  to  the  energy  gap  plus  the  mean  energy  relaxation  of  an  electron 
as  It  travels  across  the  bridge.  Now,  the  measured  dc  voltage  also 
includes  the  voltage  drops  in  the  electrodes  due  to  the  electric  field 
in  the  superconducting  electrodes.  Thus  the  peak  in  dV/dl  should  occur 
when  the  measured  dc  voltage  Is  equal  to  the  energy  gap  of  the  super¬ 
conducting  electrode  plus  the  voltage  drops  in  the  electrodes  plus  the 
mean  energy  relaxation  in  the  bridge. 

It  should  be  pointed  out  that  simple  joule  heating  In  the  sample 
would  also  raise  the  mean  energy,  thus  Increasing  the  extra  resistance. 
But  increased  heating  would  not  cause  a  decrease  In  the  occupation  of 
an  energy  level,  and  would  therefore  not  produce  a  drop  In  the  resistance 
as  the  voltage  Increased.  The  drop  in  resistance  is  a  consequence  of  the 
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length  of  the  bridge  being  comparable  to  the  Inelastic  diffusion  length. 

The  observed  peak  could  not  be  the  result  of  a  hot  spot  in  the  bridge, 
since  a  hot  spot  requires  that  as  the  temperature  is  lowered  the  peak 
should  occur  at  higher  voltages  due  to  the  need  for  increased  power  dissi¬ 
pation  to  sustain  the  hot  spot;  our  peak  occurs  at  a  temperature- 
independent  voltage  at  low  temperatures. 

The  above  analysis  indicates  that  the  width  of  the  resistance  peak 
is  a  measure  of  the  width  of  the  electron  energy  distribution  in  the 
normal  metal.  In  addition,  the  various  bumps  on  the  dV/dl  curve  should 
also  represent  structure  in  the  density  of  states  near  the  interface. 
Several  authors  3,0  ”0,3  have  calculated  that  the  occupied  density  of 
states  near  the  fermi  level  is  broadened  and,  in  addition,  the  density 
of  states  actually  oscillates  with  energy.  Thus  it  appears  that  our 
experimental  observations  are  somewhat  in  agreement  with  theory,  but  the 
agreement  is  at  best  qualitative.  A  major  difficulty  in  applying  the 
theories  is  that  all  of  the  theories  are  done  in  the  clean  limit,  while 
our  samples  are  In  the  dirty  limit. 

Rowell  and  McMillan®®  and  Rowell®®  used  a  tunnelling  geometry  to 
measure  the  density  of  states  in  clean  SN  systems,  and  they  found 
excellent  agreement  with  theory.  But  since  their  experiments  were 
tunnelling  experiments,  they  put  only  a  very  small  voltage  across  the 
normal  metal.  Our  result  appears  to  be  the  first  observation  of  oscilla¬ 
tions  In  the  density  of  states  in  dirty  systems,  as  well  as  the  first 
observation  of  the  oscillations  by  putting  a  relatively  large  voltage 
across  the  sample. 


Chapter  5 
OTHER  SAMPLES 


In  addition  to  the  lead-copper  SNS  microbridges,  other  samples  were 
produced  in  the  course  of  this  work.  Of  these  other  samples,  only  those 
with  niobium  banks  were  of  sufficient  quality  to  justify  detailed  measure¬ 
ment  of  their  electrical  properties.  But  even  these  niobium-banked  micro¬ 
bridges  performed  poorly  compared  to  the  lead-banked  microbridges.  (This 
poor  performance  was  due  primarily  to  the  relatively  long  reduced  lengths 
of  these  samples.  Because  of  the  difficulties  in  sample  preparation,  the 
niobium-banked  structures  had  relatively  thin  bridges,  of  the  order  of 
20  nm  thickness.  Thus,  in  spite  of  the  shortness  of  the  bridges  on  a 
geometrical  scale,  as  short  as  200  nm,  the  bridges  were  all  fairly  long 
in  terms  of  the  coherence  length.)  As  a  result  of  this  poor  performance, 
it  was  not  possible  to  quantitatively  compare  the  results  of  measurements 
on  these  niobium-banked  structures  to  the  available  theories  of  super¬ 
conducting  weak  links.  However,  it  was  useful  to  compare  the  electrical 
properties  of  these  niobium-banked  microbridges  to  the  results  obtained 
for  the  lead-banked  microbridges.  Thus,  in  order  to  make  this  comparison, 
the  results  on  these  other  samples  will  be  briefly  presented. 

5.1  Nb-Cu-Nb  Samples 

The  microbridges  consisting  of  niobium  banks  and  copper  bridges, 
like  the  lead-copper  microbridges,  were  SNS  microbridges.  The  current- 
phase  relation  of  these  Nb-Cu-Nb  samples  was  found  to  be  nearly  sinusoidal 
at  all  temperatures  for  which  the  CPR  could  be  measured.  Furthermore,  as 
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5.1  Log  ICR  versus  T/Tp  for  a  Nb-Cu-Nb  microbridge  (Cu  1-4). 
This  sample  had  a  bridge  about  20  nm  thick,  200  nm  wide 
and  200  nm  long. 
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Figure  5.3  I-V  data  for  sample  Cu  1-4  with  15.8  GHz  radiation. 


seen  In  Figure  5.1,  the  critical  current  was  found  to  be  nearly  exponential 
In  temperature.  This  I  versus  T  data  could  be  fit  by  Usadel's  equations, 

V 

resulting  In  a  reduced  length  i  of  20.5.  This  reduced  length  was  compara¬ 
ble  to  that  of  the  longer  lead-copper  samples. 

Figure  5.2  shows  a  set  of  I-V  curves  for  the  same  sample  CCul-4)  as 
In  Figure  5.1.  If  these  curves  are  compared  to  the  I-V  curves  for  a  lead- 
copper  sample  having  a  similar  reduced  length,  the  sets  of  curves  are  found 
to  be  quite  similar.  At  low  voltages  there  is  an  inflection  point  in  the 
I-V  curve  at  about  2  microvolts;  this  voltage  corresponds  to  a  character¬ 
istic  time  of  ‘fi/2eV  *  1.6  x  10'10  seconds.  This  time  is  in  good  agreement 

2  -10 

with  the  effective  relaxation  time  of  ■  £,  irfi/8kTc  =  1.4  x  10 
seconds.  Somewhat  above  this  inflection  point  the  I-V  curves  become 
straight  lines. 

A  relaxation  time  this  long  should  cause  hysteresis  in  the  I-V  at 
low  temperatures;  but,  as  with  the  longer  Pb-Cu-Pb  samples,  the  temperature 
could  not  be  reduced  to  a  low  enough  value  for  the  hysteresis  to  be  seen. 
However,  it  can  be  seen  in  Figure  5.2  that  at  the  lowest  temperatures 
achieved  the  I-V  curve  is  discontinuous  yet  single-valued. 

Figure  5.3  shows  I-V  curves  for  Cul-4  with  15.8  GHz  applied 
microwave  radiation.  There  are  strong  constant- voltage  steps  and  there 
are  significant  subharmonic  steps.  These  constant-voltage  steps  were 
observed  to  oscillate  for  several  periods  with  increasing  applied  rf 
power.  Data  taken  at  8.3  GHz  showed  similar  constant-voltage  steps  with 
strong  subharmonic  steps.  The  presence  of  these  subharmonic  steps  is  in 
agreement  with  an  effective  relaxation  time  of  1.4  x  10"10  seconds.  Data 
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taken  at  2.0  GHz  did  not  show  subharmonic  steps. 

It  can  be  seen  In  Figure  5.3  that  with  15.8  GHz  radiation  there  is 
a  slight  enhancement  of  the  critical  current,  and  there  is  also  a  small 
Induced  hysteresis.  For  8.3  GHz  and  2.0  GHz  radiation  no  enhancement  was 
observed  at  these  low  temperatures .  This  is  in  good  agreement  with  the 
results  obtained  for  the  longer  Pb-Cu-Pb  samples  at  low  temperatures. 

Thus  it  seems  that  the  results  of  measurements  on  the  Nb-Cu-Nb  SNS 

microbridges  support  the  ideas  presented  in  Chapter  4.  Unfortunately, 

the  long  reduced  lengths  of  these  samples  prevent  any  strong  quantitative 

comparisons  to  be  made.  It  was  not  possible  to  lower  the  temperature 

enough  to  observe  hysteresis  in  the  I-V  curves,  and  the  low  values  of 

I  made  it  impossible  to  measure  the  effects  of  applied  rf  radiation  at 
c 

any  but  the  lowest  temperatures. 

The  Nb-Cu-Nb  samples  did  have  one  strong  advantage  over  the 
Pb-Cu-Pb  samples,  however.  The  high  melting  temperatures  of  niobium  and 
copper  meant  that  prolonged  storage  of  these  samples  at  room  temperature 
did  not  alter  their  electrical  characteristics. 

5.2  Nb-Al-Nb  Samples 

Microbridges  with  niobium  banks  and  aluminum  bridges  were  produced 
in  which  the  bridge  was  about  20  ran  thick  and  one  micron  long.  Figure  5.4 
shows  a  set  of  I-V  curves  for  one  of  these  samples  (A  14-3).  These  curves 
show  no  hysteresis  or  other  noteworthy  features.  Figure  5.5  shows  the 
critical  current  versus  temperature  data  for  this  sample;  I  is  seen  to 
be  nearly  exponential  in  temperature.  This  IC(T)  relation  is  similar 
to  that  found  for  SNS  structures.  However,  the  solution  to  Usadel's 
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Figure  5.5  Log  I  versus  T  data  for  sample  A1  4-3. 
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equations  found  in  Chapter  4  is  not  applicable  to  these  samples,  since 
aluminum  is  a  superconductor  with  a  T  of  about  1  K;  there  has  not  been 
any  theory  advanced  which  predicts  an  I  versus  T  relation  for  micro- 
bridges  with  such  a  large  difference  between  the  transition  temperatures 
of  the  bridge  and  banks,  for  data  taken  near  the  transition  temperature 
of  the  bridge. 

The  current-phase  relation  for  sample  A 14-3  was  found  to  be  nearly 
sinusoidal  at  all  temperatures  of  measurement. 

Figure  5.6  shows  I-V  curves  with  2.0  SHz  applied  microwave  radia¬ 
tion.  There  are  no  subharmonic  steps,  and  there  Is  no  enhanced  critical 
current.  The  constant-voltage  steps  showed  only  a  single  oscillation 
period  with  increasing  microwave  power.  Measurements  were  also  performed 
at  7  GHz,  and  no  enhancement  of  I  was  observed. 

The  poor  performance  of  these  samples,  both  with  respect  to  the  low 
values  of  I  R  and  the  lack  of  good  constant-voltage  steps,  was  probably 
due  to  the  thinness  of  the  bridges. 

5.3  Nb-Nb-Nb  Samples 

Variable-thickness  microbridges  were  produced  in  which  both  the  bridge 
and  the  banks  were  made  from  niobium;  the  bridges  were  about  25  nm  thick 
and  the  banks  were  about  60  nm  thick. 

Figure  5.7  shows  critical -current  versus  temperature  data  for  a 
sample  (Nb  36-1 )  which  had  a  bridge  about  0.5  micron  long;  Figure  5.8 
shows  similar  data  for  a  sample  (Nb  37-2)  which  had  a  bridge  about  1.5 
microns  long.  For  both  samples  Ic  was  found  to  be  nearly  proportional  to 
(T  -T)^,  In  agreement  with  the  predictions  for  a  two-dimensional 

w 


V(^V) 


Te-T  (K) 

Figure  5.7  Log  Ic  versus  (Tc-T)  data  for  a  Nb-Nb-Nb  microbridge 
(Nb  36-1).  This  sample  had  a  bridge  about  0.5  microns 
wide  and  0.5  microns  long. 
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0.1  0.01  0.001 

I  -  T/Tc 

Figure  5.8  Log  Ic  versus  (1-T/TC)  data  for  Nb-Nb-Nb  microbridge 
(Nb  37-2) .  This  sample  had  a  bridge  about  0.5  microns 
wide  and  1.5  microns  long. 
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superconducting  film;  this  data  indicates  that  these  samples  were  extremely 

non-ideal.  Current-phase  measurements  on  Nb 36-1  found  that  the  CPR  was 

non-sinusoidal ,  becoming  multi-valued  slightly  below  Tc. 

Figure  5.9  shows  I-V  curves  for  Nb  36-1 .  There  are  two  interesting 

features.  First,  the  I-V  has  an  upward  curvature  almost  everywhere;  this 

is  probably  due  to  heating.  Second,  at  lower  temperatures  the  I-V  curve 

is  hysteretic.  Figure  5.10  shows  the  dependence  of  the  return  current 

IR  on  the  critical  current;  it  was  found  that  IR  «  Ic  .  This  indicates 

27 

that  the  hysteresis  was  due  to  heating  ,  in  contrast  to  the  relaxation¬ 
time-induced  hysteresis  seen  in  the  Pb-Cu-Pb  samples. 

These  samples  were  studied  with  applied  microwave  radiation  for 
frequencies  up  to  10  GHz,  and  only  extremely  weak  constant-voltage  steps 
were  seen.  No  enhancement  of  the  critical  current  was  observed. 


2 
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Figure  5. 
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9  I-V  data  for  sample  Nb  36-1 


Chapter  6 

SUMMARY  AND  CONCLUSION 

The  results  reported  In  this  work  represent  the  first  measurements 
of  the  electrical  properties  of  short,  high  resistance  (R  ^  1  ohm)  SNS 
microbridges.  These  small  (length  and  width  as  small  as  200  nm),  planar 
structures  were  produced  using  the  advanced  microfabrication  techniques 
of  electron-beam  lithography  and  ion-beam  etching. 

The  small  cross  section  of  the  bridges  resulted  in  a  uniform  current 
distribution  in  the  bridges  at  all  temperatures.  This  enabled  us  to  make 
the  first  experimental  determination  of  the  critical  current  versus  tempera 
ture  relation  at  temperatures  far  from  the  critical  temperature  Tc  of  the 
banks;  the  Ic  versus  T  relation  was  measured  for  all  temperatures  between 
T  and  0.05  T/T  .  We  found  that  our  experimental  I  (T)  relation  was  in 

V  C 
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good  agreement  with  theoretical  predictions  based  on  Usadel's  equations. 
(Usadel's  equations  are  a  reformulation  of  Gorkov's  formulation  of  super¬ 
conductivity.  They  are  valid  for  any  superconducting  system  in  the  dirty 
limit,  not  just  for  SNS  microbridges.  This  work  has  provided  the  first 
experimental  verification  of  these  equations.)  This  successful  fit  of 
theory  to  experiment  produced  a  reduced  length  l  for  each  sample,  with 
l  *  L/ (fiv^A/ 6irkTc)is  where  L  is  the  physical  length  of  the  bridge  and  A 
Is  the  electron  mean  free  path  in  the  bridge.  Our  samples  had  reduced 
lengths  ranging  from  8.9  to  30. 

The  low-voltage  (V  <  TO  I  R)  current-voltage  (I-V)  characteristics 

1# 

were  studied  in  detail,  and  we  observed  for  the  first  time  the  effects  of 
an  intrinsic  relaxation  time  on  the  I-V  curves.  We  found  that  these  I-V 
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curves  could  be  characterized  by  an  Intrinsic  temperature-independent 
relaxation  time  »  l  irfi/8kTc  where  l  is  the  reduced  length  of  the 
bridge  determined  from  the  IC(T)  data.  We  found  that  at  high  tempera¬ 
tures  there  was  an  inflection  point  in  the  I-V  curve  at  a  voltage 
V  ,'''h/2eTe^.  At  low  temperatures  the  I-V  curves  were  hysteretic,  with 
the  I-V  curve  returning  to  the  zero-voltage  axis  at  a  return  current  IR. 

At  all  temperatures  the  I-V  curve  became  a  straight  line  at  a  voltage 

u 

lower* 

We  observed  that  at  low  temperatures  IRR  and  VlQwer  were  independent 

of  temperature,  with  both  IRR  and  V^ower  having  the  same  low- temperature 

value  of  ^rR/2exeff  (rR  ^5).  We  found  that  rR  was  nearly  independent  of 

the  sample  length.  These  observations  confirmed  both  the  temperature 

2 

independence  of  xfiff  and  the  %  dependence  of  teff. 

This  was  the  first  observation  of  hysteresis  in  SNS  structures.  It 
was  also  the  first  time  that  the  hysteresis  in  any  type  of  planar  micro¬ 
bridge  could  be  definitely  attributed  to  a  relaxation  time  rather  than 
heating  or  capacitance. 

It  was  found  that  the  nature  of  the  I-V  curve  with  applied  microwave 
radiation  was  also  dependent  on  the  effective  relaxation  time.  These 
I-V  curves  showed,  for  some  microwave  frequencies,  a  microwave-enhanced 
critical  current.  It  was  determined  that  if  the  microwave  period  rrf 
was  shorter  than  either  fi/2eIcR  or  t^,  the  critical  current  was  enhanced. 
If  tr^  was  longer  than  both  of  these  times,  there  was  no  enhancement. 

This  result  indicates  that  the  popular  gap-enhancement  explanation  for  a 
microwave-enhanced  critical  current  does  not  apply  to  these  SNS 
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microbridges,  but  that  the  enhancement  is  more  likely  caused  by  a 
microwave-induced  spatial  broadening  of  the  quasiparticle  distribution. 

Constant-voltage  steps  were  observed  in  the  I-V  curves  for  all 
frequencies  of  applied  microwave  radiation  used.  Subharmonic  steps  were 
observed  whenever  the  criterion  for  an  enhanced  critical  current  was  in 
effect,  l.e. ,  for  less  than  or  fi/2eIcR.  If  was  longer  than 
both  of  these  times,  there  were  no  subharmonic  steps.  The  good  thermal 
properties  of  these  microbridges  was  demonstrated  by  the  observation  of 
constant-voltage  steps  at  dc  voltages  In  excess  of  2  mV  for  the  shortest 
samples. 

The  high  resistance  of  these  samples  allowed  us  to  make  the  first 
measurements  of  the  I-V  curves  for  voltages  much  greater  than  the  energy 
gap  of  the  electrodes.  A  high-voltage  peak  in  the  dV/dl  curves  was 
observed,  and  it  was  speculated  that  this  peak  is  due  to  an  extra  dissi¬ 
pation  in  the  electrodes.  In  addition,  the  low-voltage  resistance  was 
found  to  have  an  extra,  temperature-dependent  component.  The  temperature 

dependence  of  this  resistance  was  found  to  agree  with  the  phenomenologi- 

50 

cal  theory  of  Pippard  et  ajk  ,  a  theory  which  attributes  the  extra 
resistance  to  energy  dissipation  in  the  superconducting  electrodes. 


APPENDIX 


A. I  Numerical  Solution  of  Usadel's  Equations  for  SNS  Structures 

It  was  pointed  out  in  Chapter  4  that  the  correlation  functions  for 
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an  SNS  structure  obey  Usadel's  equations: 

•  o  fc  fd  -  ft"  &  *  \  -1 V*  TT 

v.  u  "  n  ' 

where  D  *  -j  vfA  and  w  *  (2n+l )irkT/fi.  Note  that  there  are  an  Infinite 

number  of  F  's. 
n 

The  current  through  the  bridge  is  given  by 

M*3!1  I  V 

s  e  n>0  n 


where 


m 


|SH(TJ 


“  7  rn  3>T  n  ~dx 


Likharev'1  has  demonstrated  that  these  equations  can  be  solved  in 
the  following  manner.  First,  use  the  coherence  length  as  a  scale  of 
spatial  dimension,  i.e.,  units  of  ('fiv^A/eirkT)  .  Second,  write  the 
complex  function  F  as  F  a  sin0e1Y  where  9  and  y  are  real. 

Let  8q  be  the  value  of  0  at  the  center  of  the  bridge,  and  0B  the 
value  of  0  at  the  banks.  Then  we  have  the  following  equations: 

4|  *  In  f — 7 - V]  +  2(2n+l)(cos9  -  cose)  *  , 

"  n  Uin20c  sinV  0 


1  *  cHT  *  2  d073x  * 
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<<>  a  phase  difference  across  bridge  *  2j 


where 


l»0S  s1"  9 


0B  a  tan"^  (2n+T~{lkT)  A(T)  the  value  °f  the  energy 

gap  In  the  superconducting  electrodes. 

The  difficulty  with  the  above  equations  Is  that  we  are  Interested  In 
knowing  j  as  a  function  of  2  and  <j>.  These  integral  equations  were  solved 
using  the  following  iterative  method. 

First,  by  trial  and  error,  values  of  0Q  and  jn  were  found  which  gave 
the  desired  2  at  $  *  0.05.  (With  a  little  experience  this  was  fairly 
easy.  The  integrals  were  evaluated  using  Simpson's  rule  with  a  variable 
spacing;  the  spacing  between  points  was  varied  as  1/9  ,  and  usually  1000 
points  were  used.) 

Once  these  initial  values  were  determined,  the  full  jn  versus  <p 
relation  was  found  by  incrementing  <J>  by  0.05  and  finding  the  values  of 
jn  and  0Q  which  gave  the  desired  2.  (This  procedure  was  made  automatic 
by  using  Newton's  method  for  finding  the  zeros  of  a  function.) 

This  procedure  gave  the  current  jn  as  a  function  of  <p.  For  all  the 
lengths  and  temperatures  applicable  to  our  samples,  jn  was  found  to  vary 
as  sin$  within  10%. 

The  above  procedure  had  to  be  repeated  for  increasing  values  of  n 


until  increasing  n  changed  Zjn  by  less  than  .01%.  (For  T/Tc  »  0.6  only 
n  *  1  was  needed;  for  T/T  *  0.1  values  of  n  up  to  4  were  required.) 
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A. 2  Numerical  Simulation  of  the  Shunted  Junction  Model 
The  reslstlvely  shunted  junction  model  is  given  by 

!dc  +  *rf  s1n^rft^  a  s1n<{>  +  V/R 
and 

V  -  ft/2e)  ||  . 

If  currents  are  measured  In  units  of  I  ,  voltage  in  units  of  I  R,  and 

C  V 

time  in  units  of  fi/2eIcR,  then  these  equations  become 

*dc  +  *rf  sin(u»rft)  *  si n<|>  +  V 
and 


These  equations  were  solved  numerically  by  stepwise  Integration. 

That  Is,  initially  V  *  0,  <J>  3  0,  and  t  *  0.  The  time  was  then  incre¬ 

mented  by  At,  and  I$,  V,  and  <j>  were  recalculated  by 

I$  3  s1n<j> 

V  =  Xdc  +  lrf  s1nKft)  *  h 
<j)  »  VAt  . 

This  process  was  continued  for  several  periods,  with  the  values  of  I$, 

V,  and  <j)  recorded  at  each  step.  The  end  result  was  Is(t),  V(t),  and 

<|>(t)  data.  Typically,  At  ■  .01,  and  the  values  of  I$,  V,  and  <(»  were 

not  recorded  until  several  periods  had  passed  (about  1000  steps)  to 
ensure  that  the  Iteration  had  stabilized. 

In  order  to  calculate  an  I-V  curve,  the  above  procedure  was  carried 
out  for  several  values  of  !dc'  Vdc  was  determined  by  averaging  V  over 
a  few  hundred  periods. 
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The  effects  of  a  characteristic  relaxation  time  rre-j 
Incorporated  by  calculating  I$  with  the  equation 


ax 


2  1 


I$  *  s1n<j> 


rel  ax 


where 


^ relax  "  ^^e^c^Trelax^  ' 


were 
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